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SUMMARY 

In  this  report  we  show  how  proper  cooperation  among  radio  nodes  may  provide  diversity  gain,  also  for 
single  antenna  systems.  We  consider  first  the  connectivity  of  a  wireless  network  and  show  how  it  can 
henefit  from  cooperation.  Then,  we  consider  some  specific  forms  of  cooperations,  based  on  distributed 
space-time  coding,  in  both  single  and  multi-user  contexts.  Finally,  we  pay  a  special  attention  to  the  case 
where  the  final  destination  has  a  multi-antenna  receiver.  In  such  a  case,  we  may  establish  a  virtual  MIMO 
link  between  the  relays  and  the  final  destination,  which  makes  possible  to  benefit  also  from  the  MIMO 
spatial  multiplexing  gain. 
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One  of  the  distinguishing  features  of  radio  channels  is  its  fading  nature.  Channel  fading  has  received 
a  considerable  attention  as  it  can  cause  severe  performance  losses.  Since  fading  is  a  multiplicative 
phenomenon,  the  best  way  to  counteract  it  is  not  a  simple  increase  of  transmission  power.  The  other 
major  characteristic  of  radio  propagation  is  its  broadcasting  nature:  The  information  emitted  by  a  radio 
source  travels  in  fact  in  all  directions,  even  though  with  possibly  different  amplitudes,  depending  on 
the  antenna  directivity.  This  means  that  part  of  the  transmitted  energy,  a  typically  scarce  resource,  is 
inevitably  wasted.  To  limit  such  a  waste,  one  should  use  directive  antennas,  but  sometimes  this  could  be 
highly  impractical,  especially  on  portable  handsets,  as  it  requires  the  use  of  large  antennas  (with  respect 
to  the  transmit  wavelength). 

The  most  effective  strategy  to  combat  channel  fading  is  diversity.  Diversity  relies  on  the  existence  of  more 
than  one  path  between  transmitter  and  receiver  through  which  the  information  bits  may  be  conveyed.  If 
the  paths  are  random,  but  not  strongly  correlated  to  each  other,  one  can  find  optimal  strategies  to  combine 
the  multiple  replicas  of  the  received  signal  in  order  to  minimize  the  bit  error  rate  (BER).  The  most  known 
form  of  diversity  is  spatial  receive  diversity,  obtainable  when  the  receiver  has  multiple  receive  antennas. 
Nevertheless  the  last  years  have  witnessed  a  huge  amount  of  research  on  methods  capable  of  achieving  a 
spatial  diversity  even  in  the  case  where  the  receiver  has  only  one  antenna,  but  the  transmitter  has  multiple 
antennas.  In  such  a  case,  it  is  necessary  to  transmit  the  data  using  a  particular  form  of  coding  that  spreads 
the  information  bits  across  the  two-dimensional  space-time  domain,  using  the  so  called  space-time  coding 
[6],  [7].  If  properly  designed,  space-time  coding  guarantees  full  diversity  gain,  even  in  cases  where  the 
transmitter  does  not  have  any  information  about  the  channels. 

However,  if  both  transmitters  and  receivers  have  only  one  antenna,  it  looks  like  there  should  be  no  way 
to  achieve  any  spatial  diversity.  Indeed,  this  is  not  true,  if  more  radio  nodes  (relays)  contribute  to  the 
transmission  from  one  source  to  the  intended  destination.  In  fact,  it  is  precisely  the  broadcasting  nature 
of  radio  transmissions  that  makes  possible  the  propagation  of  the  same  information  through  several  relay 
nodes  willing  to  cooperate  with  the  original  source.  If  the  radio  nodes  that  have  correctly  decoded  the 
same  message  cooperate  with  each  other  to  re-transmit  the  information  to  the  final  destination,  the  overall 
system  may  benefit  of  spatial  diversity,  even  if  each  radio  node  has  only  one  antenna.  This  technical 
report  focuses  on  this  idea  and  on  the  ways  to  achieve  diversity  through  cooperation. 

Multihop  radio  networking  is  a  research  field  whose  study  started  long  time  ago  (see,  e.g.  [1]  and  the 
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references  therein).  Some  basic  theorems  on  the  capacity  of  relaying  networks  were  established  in  [2],  [8], 
[14],  [15],  [4],  [5],  [3].  Conventional  relaying  can  indeed  be  seen  as  a  particular  form  of  distributed  space- 
time  coding,  where  the  same  information  is  transmitted  from  different  (spatial)  locations,  at  different  times. 
From  this  perspective,  relaying  can  thus  be  interpreted  as  a  repetition  code  in  the  space-time  domain.  As 
well  known,  repetition  coding  is  not  the  best  coding  strategy  and  thus  one  should  achieve  a  considerable 
gain  by  using  more  sophisticated  space-time  coding  techniques.  This  form  of  space-time  coding  that 
coordinates  the  transmission  of  source  and  relays  is  called  distributed  space-time  coding.  Quite  recently, 
the  interest  about  relaying  networks,  especially  in  the  form  of  cooperation  among  nodes,  has  increased 
considerably.  One  result  that  sparkled  great  interest  was  that  cooperation  among  users  can  increase  the 
capacity  in  an  uplink  multiuser  channel  [16].  A  thorough  analysis  of  the  diversity  gain  achievable  with 
cooperation  was  given  in  [17],  [18],  [19],  where  different  distributed  cooperation  protocols  were  compared. 
Cooperation  was  proved  to  be  very  useful  to  combat  shadowing  effects,  as  shown  in  [20],  and  it  can 
occur  in  different  forms,  as  suggested  in  many  recent  works,  like  e.g.  [18],  [21],  [22],  [23],  [24],  [25], 
[26],  [27],  [28],  [29],  [30].  Among  all  these  possibilities,  there  is  the  basic  idea  of  using  the  relays  as 
if  they  were  the  antennas  of  a  multi-antenna  transmitter.  Within  this  perspective,  cooperation  induces  a 
virtual  array  between  transmitter  and  receiver.  This  opens  the  field  to  a  huge  amount  of  possibilities, 
provided  by  all  coding  methods  devised  for  multi-input/multi-output  (MIMO)  transceivers.  In  this  more 
general  framework,  the  relaying  problem  becomes  the  problem  of  mapping  the  algorithms  valid  for  real 
MIMO  systems  into  cooperation  strategies  of  a  virtual  MIMO  system. 

This  technical  report  is  organized  as  follows.  In  Section  ??,  we  study  the  connectivity  of  a  wireless 
network  and  show  how  a  single  user  system  can  achieve  a  considerable  gain  thanks  to  the  cooperation  with 
relay  nodes  scattered  randomly  in  a  given  territory.  In  Section  II,  we  introduce  the  so  called  distributed 
space-time  coding  (DSTC)  strategy,  where  source  and  relays  transmit  in  a  coordinate  manner  according  to 
a  space-time  code  distributed  among  the  cooperating  nodes.  In  Section  III,  we  consider  in  detail  a  specific 
single  user  system  where  source  and  relay  adopt  a  block  distributed  Alamouti  strategy,  for  transmission 
over  frequency  selective  channels.  We  show  the  performance  achieved  over  simulated  as  well  as  real 
data.  Although  DSTC  follows  the  main  ideas  of  conventional  space-time  coding  (STC)  techniques,  there 
are  aspects  that  clearly  differentiate  DSTC  from  STC.  The  main  differences  are  related  to  the  fact  that  in 
DSTC  the  cooperating  antennas  are  not  co-located.  This  implies  that  the  signals  received  by  the  destination 
from  source  and  relays  might  arrive  out  of  synchronization  and  without  the  right  power  balance  between 
them.  We  analyze  these  aspects  in  Section  III-C,  where  we  consider  possible  ways  to  distribute  the  power 
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between  source  and  relays  optimally,  in  order  to  minimize  the  average  BER,  and  in  Section  III-F,  where 
we  consider  the  synchronization  problem. 

The  performance  of  a  cooperative  network  improves  as  the  density  of  the  relay  nodes  increases.  However, 
taking  into  account  both  technical  as  well  as  economic  factors,  this  implies  that  a  cooperative  strategy 
is  appealing  only  if  the  cost  of  the  relay  nodes  is  extremely  small.  One  possibility  is  to  use  as  relays, 
in  cellular  networks,  portable  phones  that  are  in  standby,  whose  owners  have  previously  agreed  on 
their  availability  to  act  as  relays,  maybe  against  appropriate  incentives.  The  other  possibility  consists  in 
deploying  extremely  simple  relays,  such  as  for  example  non-regenerative,  or  amplify-and-forward  (A&F), 
relays  that  amplify  and  forward  the  received  messages.  Such  relays  need  only  the  radio-frequency  section 
(antenna  and  amplifier)  and  are  then  much  more  economical  than  any  other  transceiver.  To  assess  the 
performance  of  such  relays,  in  Section  IV  we  compare  the  performance  of  regenerative,  i.e.  decode- 
and-forward  (D&F),  relays  and  non-regenerative  relays.  In  Section  V  we  compare  alternative  DSTC 
strategies  for  a  multi-user  system.  Finally,  in  Section  VI,  we  derive  the  outage  probability  of  a  relay 
scheme  composed  of  the  cascade  of  a  broadcasting  channel  followed  by  a  virtual  MIMO  channel. 

I.  On  the  connectivity  of  cooperative  vs.  non-cooperative  wireless  communication 

NETWORKS 

In  a  seminal  work,  Gupta  and  Kumar  derived  the  conditions  for  the  asymptotic  connectivity  of  a 
network  composed  of  nodes  uniformly  distributed  over  a  disc  of  unit  area,  as  the  number  of  nodes  goes 
to  infinity  [8].  In  this  section,  we  incorporate  the  channel  fading  and  we  provide  the  conditions  for 
the  network  connectivity,  in  case  of  single  or  multi-antenna  transceivers.  In  particular,  we  derive  closed 
form,  albeit  approximate,  expressions  for  the  spatial  density  with  which  the  nodes  must  be  deployed  in 
order  to  insure  the  network  connectivity  with  a  desired  probability.  Finally,  we  show  how  to  improve  the 
connectivity  by  allowing  nearby  nodes  to  transmit  in  a  cooperative  manner,  using  a  distributed  space-time 
coding  strategy,  in  order  to  get  spatial  diversity  gain. 

One  of  the  fundamental  issues  in  the  design  of  a  wireless  network  is  the  minimal  transmit  power  that 
guarantees  the  network  connectivity,  i.e.  the  property  that  each  node  is  connected  to  each  other  node.  In 
a  seminal  work,  Gilbert  modeled  the  network  topology  as  a  homogeneous  two-dimensional  Poisson  point 
process,  with  constant  density  A,  where  two  points  are  linked  if  their  distance  is  less  than  a  coverage 
radius  r(n)  that  depends,  on  its  turn,  on  the  transmit  power  as  well  as  on  the  number  of  nodes  n  [9]. 
Starting  from  this  simple,  yet  meaningful,  model  and  applying,  for  the  first  time,  basic  concepts  from 
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continuum  percolation,  Gilbert  showed  that  there  exists  a  critical  value  Ac  of  the  node  density  such  that, 
for  A  <  Ac  the  network  is  disconnected  whereas  for  A  >  Ac  there  is  a  nonnull  probability  that  the  network 
contains  a  giant  component  composed  of  an  infinite  number  of  nodes,  thus  creating  the  possibility  for 
long  distance  communications  through  multiple  hops  [9].  Of  course,  the  existence  of  a  component  with 
infinite  nodes  does  not  guarantee  the  network  connectivity,  as  some  nodes  could  be  isolated  from  that 
giant  component.  The  connection  between  continuum  percolation,  covering  strategies  and  the  geometry 
of  wireless  networks  has  been  further  investigated  quite  recently  in  [10]. 

In  another  fundamental  paper,  Gupta  and  Kumar  derived  the  minimum  (critical)  power  that  guarantees 
the  network  connectivity  with  probability  one  [8].  More  specifically,  Gupta  and  Kumar  proved  that  if  the 
network  is  composed  of  a  set  of  n  nodes  randomly  distributed  over  a  disc  of  unit  area  and  each  node 
transmits  with  a  power  that  allows  the  coverage  of  a  circle  of  radius  r(n)  that  scales  with  the  number 
of  nodes  n  as  follows 


the  network  is  asymptotically  connected  with  probability  one  if  and  only  if  c(n)  tends  to  infinity  as  n 
goes  to  infinity. 

In  this  work  we  assume,  as  in  Gilbert’s  work,  a  homogeneous  Poisson  point  process  and  we  find  the 
coverage  radius  ro(n)  that  guarantees  that  there  are  no  isolated  nodes,  where  the  connection  between  nodes 
is  defined  in  terms  of  out-of-service  probability,  with  respect  to  a  target  BER,  for  a  given  channel  fading 
model.  Then,  exploiting  Penrose’s  theorem  [48]  on  the  relationship  between  the  connectivity  and  the 
degree  of  a  random  geometric  graph,  we  derive  asymptotic  results  on  the  network  connectivity.  Building 
on  these  derivations,  we  show  how  the  connectivity  improves  by  using  multiantenna  terminals.  Finally, 
we  extend  the  analysis  to  cooperative  networks  where  nearby  nodes,  i.e.  nodes  within  the  coverage  of 
each  other,  transmit  in  a  coordinated  manner  according  to  a  distributed  space-time  coding  strategy,  in 
order  to  achieve  transmit  diversity  gain. 

A.  Connectivity  of  a  random  geometric  graph 

We  start  reviewing  some  results  on  the  connectivity  of  geometric  random  graphs,  as  derived  by 
Bettstetter  in  [47].  We  assume  that  the  nodes  are  distributed  according  to  a  two  dimensional  (2D)  Poisson 
point  process,  with  constant  spatial  density  p.  One  of  the  key  properties  of  this  kind  of  point  process 
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is  that  the  pdf  of  the  distance  between  a  point  and  its  nearest  neighbor  is  a  Rayleigh  pdf.  In  formulas, 
denoting  with  H  the  random  variable  indicating  the  distance  between  one  node  and  its  nearest  neighbor, 
the  pdf  of  S  is 

Pe{0  =  (2) 

where  p  is  density  of  the  Poisson  process.  Starting  from  (2),  and  denoting  by  ro(n)  the  coverage  radius 
of  each  node,  the  probability  that  a  node  is  isolated  is 

^isolated  =  1  -  <  ro(n)}  =  (3^ 

The  probability  that  there  are  no  isolated  nodes  or  that,  equivalently,  the  degree  of  the  graph  is  strictly 
greater  than  zero,  is,  approximately  [47] 

(4) 

This  expression  is  approximated  as  it  implicitly  assumes  that  the  events  for  which  the  nodes  are  not 
isolated  are  statistically  independent  of  each  other,  which  is  not  true,  in  general.  Nevertheless,  the  validity 
of  this  approximation  was  tested  in  [47]  and  we  will  also  test  it  through  simulation  results  later  on.  Indeed, 
the  validity  of  the  approximation  depends  on  the  product  pr^ln).  If  prQ{n)  <C  1,  the  approximation  is 
clearly  not  valid  at  all.  Conversely,  if  prQ(n)  1,  the  approximation  is  more  and  more  valid.  On  the 
other  hand,  since  we  are  interested  in  the  event  that  the  network  is  connected  with  high  probability,  i.e. 
Peon  is  close  to  one,  this  implies  that  we  are  interested  in  the  case  where  the  product  prQ{n)  is  high. 
This  means  that  we  can  use  the  above  approximation  in  cases  of  practical  interest.  Within  the  validity 
of  the  approximation  leading  to  (4),  inverting  (4),  the  coverage  radius  vq  that  ensures  that  the  degree  of 
the  network  is  greater  than  zero  is: 


Of  course,  this  radius  does  not  guarantee  the  network  connectivity,  as  the  network  could  be  composed 
of  isolated  clusters,  with  no  isolated  nodes.  In  general,  in  fact,  the  radius  that  guarantees  the  global 
connectivity  is  greater  than  the  radius  that  simply  guarantees  that  the  degree  of  the  network  is  strictly 
greater  than  zero.  Nevertheless,  Penrose  proved  that  there  is  a  basic  relationship  between  connectivity 
and  degree  of  a  geometric  graph  [48].  In  particular,  Penrose  proved  that,  denoting  with  r{K  >  k)  the 
minimum  distance  tq  such  that  a  random  graph  is  fe-connected  and  with  r{6  >  k)  the  minimum  distance 
such  that  the  graph  has  minimum  degree  k,  then  [48] 

lim  Pr{r(K  >  k)  =  r{5  >  k)}  =  1,  VA:.  (6) 

n— >oo 
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This  means  that,  asymptotically  as  n  tends  to  infinity,  the  coverage  radius  that  insures  a  degree  k 
tends  to  coincide  with  the  radius  that  yields  connectivity  k.  Hence,  the  equivalence  (6)  guarantees  that, 
asymptotically,  as  n  tends  to  infinity,  the  radius  ro(n)  in  (5)  tends  to  coincide  with  the  radius  that 
guarantees  also  the  connectivity.  Hence,  the  value  ro(n)  in  (5)  denotes,  approximately,  the  minimum 
coverage  that  each  node  has  to  provide  in  order  to  guarantee  the  connectivity  of  the  whole  graph,  with 
a  given  prohahility. 


Fig.  1.  Probability  of  connectivity  vs.  SNRt  reported  at  the  transmit  side  (pr/No)  for  a  network  composed  of  n  =  1000 
nodes  uniformly  distributed  over  a  square  of  200  x  200m^. 


result,  recalled  in  (1).  On  the  contrary,  we  show  next  that  (5)  and  (1)  are  indeed  in  good  agreement.  In 
fact,  if  we  rewrite  the  connectivity  prohahility  Pcon{n)  as  Pcon{n)  =  1  —  e(n),  with  e(n)  <C  1,  we  may 
use  the  first  order  Taylor  series  expansion  of  pcon{n)^^'^  as  Pcon{n)^^^  1  —  e{n)jn.  As  a  consequence, 
(5)  becomes,  approximately. 

Since  in  Gupta  and  Kumar  the  node  density  is  p  =  n,  inserting  this  value  in  (7),  we  can  verify  that  (1)  is 
in  perfect  agreement  with  (5)  if  we  set  c(n)  =  —  log(e(n)).  In  fact,  requiring  asymptotic  connectivity  with 
probability  one  implies  that  lim  e(n)  =  0  and  this  corresponds  to  the  condition  lim  c{n)  =  oo  found 

n— >oo  n— >cxD 

in  [8].  Interestingly,  the  equivalence  between  these  two  different  derivations  provides  also  a  meaning  to 
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the  term  c(n)  appearing  in  [8]. 

A  test  of  validity  of  the  approximations  leading  to  (4)  is  reported  in  Fig.  1,  that  shows  pcon  as  a  function 
of  the  SNR  (referred  to  the  transmit  side),  defined  as  SNRt  :=  Pt/Nq,  as  given  in  (4)  (solid  line),  and 
the  prohahility  that  the  network  is  connected  (dashed  line),  estimated  over  200  independent  realizations 
of  random  geometric  graphs  composed  of  1000  points  distributed  over  a  toroidal  surface,  to  avoid  border 
effects.  We  can  see  a  very  good  agreement  between  theory  and  simulation.  The  slight  difference  is  due 
to  the  fact  that  the  independence  assumption  is  not  exactly  valid  and  that  (6)  is  valid  only  for  n  going  to 
infinity.  It  is  interesting  to  notice,  from  Figure  1,  the  rapid  change  from  probability  zero  to  probability 
one.  This  is  indeed  a  characteristic  of  random  geometric  graphs,  that  is  reminiscent  of  phase  transitions 
in  chemistry. 


If  we  wish  to  increase  the  fault  tolerance  of  the  graph,  we  need  to  increase  the  connectivity  order.  It  was 
shown  in  [47]  that  the  probability  that  a  graph  is  fe-connected  is  approximately  equal  to  the  probability 
that  the  minimum  degree  of  the  graph  is  at  least  A:;  that  is 


Prfc  :=  Pr{G  is  k  —  connected}  ~  Fi{dmin  >  k} 


k-l 


-pnr^ 


i=0 


(8) 


This  expression  does  not  admit  an  inverse  in  closed  form,  but  it  is  certainly  invertible  as  it  is  a  monotonic 
increasing  function  of  vq.  The  value  of  tq  providing  the  desired  value  Vk,  to  be  found  numerically, 
guarantees  the  /c-connectivity,  within  the  limits  of  approximations  in  (8). 


B.  Connectivity  of  a  Wireless  Network 

In  a  wireless  network,  besides  the  position  of  the  radio  nodes,  there  is  one  more  source  of  randomness: 
the  fading  of  the  radio  links.  The  connectivity  of  a  wireless  network  is  then,  in  general,  much  more  difficult 
to  study  than  the  connectivity  of  a  random  graph.  An  interesting  recent  approach  is  proposed  in  [13] 
(see  also  the  references  therein).  In  this  section,  we  propose  an  alternative  definition  of  connectivity  of 
a  wireless  network.  In  case  of  fading,  it  is  in  fact  necessary  to  start  with  the  definition  of  connectivity. 
We  introduce  the  out-of-service  probability  Pout,  defined  as  the  probability  that  the  bit  error  rate  (BER) 
exceeds  a  target  BER,  let  us  say  Pq,  and  we  say  that  two  nodes  are  linked  to  each  other  if  the  out- 
of-service  probability  on  their  link  does  not  exceed  the  prescribed  value  Pout-  This  is  equivalent  to  say 
that  the  BER  on  the  link  may  exceed  the  target  BER  for  no  more  than  a  percentage  Pout  of  time. 
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In  the  following  sections,  we  show  how  to  choose  the  coverage  radius  in  order  to  accommodate  for 
such  a  system  requirement.  To  simplify  the  theoretical  analysis,  we  assume  that  the  channels  between 
different  pairs  of  nodes  are  independent  Rayleigh  flat-fading  channels,  with  variance  a\,  and  that  all 
links  are  interference-free.  This  means  that  there  is  an  implicit  MAC  protocol  that  guarantees  the  users’ 
separability. 


C.  Connectivity  of  a  SISO  Flat-Fading  Network 

Using  QAM  constellations,  the  bit  error  rate  on  a  SISO  flat  fading  channel  is 

Pe  =  CMQ  <  CM  (9) 

where  Sb  is  the  energy  per  bit,  Nq  is  the  noise  variance,  h  is  the  flat-fading  coefficient,  and  cm  and  qm 
are  two  coefficients  that  depend  on  the  order  M  of  the  QAM  constellation  as  follows  [32] 

^  v/M-1 

Cm  =  4  — - , 

\fM  •  log2  M 

9M  =  j^^log2M.  (10) 

We  assume  that  the  channels  are  Rayleigh  fading,  so  that  |/ip  is  an  exponential  random  variable  with 
expected  value  =  l/r“,  where  r  is  the  link  length  and  the  exponent  a  depends  on  the  environment 
where  the  propagation  takes  place.  Typically,  a  is  between  two  and  six.  Exploiting  the  upper  bound  in 
(9),  we  can  upper-bound  the  out-of-service  probability  as 

Pout=  Pr{Pe  >  Po}  <  PT{cMe-3-^^\^\P^o  >  Po}.  (11) 


Hence  Pout  can  be  upper-bounded  as 

<  1  —  e-^O^Og{cM/Po)/{9M£b(Tl) 


(12) 


We  say  that  a  link  is  reliable,  and  it  is  then  established,  if  the  out-of-service  event  occurs  with  a  probability 
smaller  than  a  given  value.  Setting  =  1/ra  in  (12)  and  inverting  (12),  we  find  the  coverage  radius 


Tcn-v  — 


gM^blog^l  -  Pout) 


I/q 


Ao  log  (cm /-Po) 

Since  Pout  is  typically  small,  we  can  use  the  approximation  log(l  —  Pout) 


(13) 

—Pout  to  rewrite  (13)  as 


Pov 


/"'w' 


dM^bPout 

Aolog(cM/Po) 


(14) 
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Equating  (14)  to  (5),  we  get  the  relationship  between  the  node  density  and  the  transmitted  power  necessary 
to  insure  the  network  connectivity,  for  a  given  number  of  nodes  n.  For  example,  the  minimum  SNR  that 
guarantees  the  network  connectivity,  with  probability  pcon,  for  a  given  node  density  p,  is: 


No 


\og{cM/Po)  f  -log(l  -pcon)\  ^ 
9  M  Pout  \  P  I 


(15) 


As  expected,  we  see  that  to  achieve  a  desired  probability  of  connectivity,  we  can  decrease  the  transmitted 
energy  if  we  increase  the  node  density  p^.  At  the  same  time,  the  transmitted  energy  must  increase  when  a 
increases.  In  a  sensor  network  scenario,  where  it  is  important  to  foresee,  a  priori,  the  density  with  which 
the  nodes  should  be  deployed,  in  a  given  area,  in  order  to  guarantee  a  certain  probability  of  connectivity, 
(15)  can  be  rewritten  by  making  explicit  the  minimum  node  density  p  that  guarantees  the  connectivity, 
for  a  given  set  of  sensors  having  a  specified  transmitted  power 

-log(l-pcon)  /No  log(cM/-Po)\^^“ 


Pmin  — 


TT 


(16) 


D.  Connectivity  of  a  MIMO  Flat-Fading  Network 

We  show  now  how  the  connectivity  improves  if  the  radio  nodes  have  multiple  antennas.  In  such  a  case, 
the  network  may  benefit  from  the  diversity  gain.  To  make  a  fair  comparison  with  the  single  antenna  case 
seen  before,  denoting  with  ut  the  number  of  transmit/receive  antennas,  we  set  the  power  transmitted  by 
each  antenna  equal  to  px/nr,  where  px  is  the  overall  transmit  power. 


Using  an  nx  x  nx  MIMO  system,  using  a  space-time  coding  technique  capable  of  achieving  full  diversity, 
the  error  probability  is 


Pe  =  CM  Q 


/ 

V 


<  CM 


(17) 


having  introduced,  in  the  last  approximation,  the  random  variable  z  :=  '}2’i=i  The  out-of-service 
probability  can  then  be  upper  bounded  as  follows 


T’out  =  Dz 


f  Nqux  log{cM/Po)\ 


QM^b 


'We  must  keep  in  mind,  however,  that  this  result  has  been  obtained  by  assuming  lack  of  interference. 
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where  Dz{z)  is  the  cumulative  distrihution  function  (CDF)  of  z.  If  the  channels  are  Rayleigh  fading, 


independent,  and  all  with  the  same  variance  ct?,  the  CDF  of  z  is 


n%  —  l 


Dz{z)  =  1  - 


fc=0 


,  k  1 

Z  \  1 


kl' 


For  small  values  of  z,  more  specifically  for  2:  <C  cr^,  Dz{z)  can  be  approximated  as 


Dz{z) 


(Jt 


'h/  ''^T- 

For  the  sake  of  finding  closed  form,  albeit  approximated,  expressions,  it  is  useful  to  introduce  the 


1 


(18) 


(19) 


normalized  random  variable  x  =  zla\,  whose  CDF  is 

Dx{x)  =  1  -  e-" 


""  k\^^V 

k=0  ^ 


(20) 


Repeating  the  same  kind  of  derivations  as  in  the  SISO  case,  the  out-of-service  probability  for  the  MIMO 
case  can  be  written  as 


-Pout  <  Dx 


/  NpUT  log(cM/Po)\ 

V  gnSbcrl  J  ' 


From  (21),  setting  =  1/ra,  we  can  derive  the  coverage  of  each  node^ 


T ccvj  — 


gM^b 


D-x\P. 


1  1/a 


out  J 


(21) 


(22) 


NqIIt  log(cM/Po) 

To  derive  an  approximate  closed  form  expression  for  since  we  are  interested  in  small  values  of  the 
out-of-service  probability,  we  can  use  the  approximation  (20)  to  invert  Dx{x).  The  result  is 


_ _ (nl\p 

[NonTlog{cM/Po)  ^  ^ 


1/q 


(23) 


Equating  (22)  to  (5),  we  get  the  minimum  transmit  power  guaranteeing  the  connectivity  of  a  MIMO 
network 


P  _  nrlog(cM/Po) 
^0  gM{n‘^'-PoutyP^^ 


-  log  (^1  -  Pcon^ 
irp 


(24) 


As  a  numerical  example,  in  Figure  2  we  show  the  density  p,  as  a  function  of  the  transmitted  energy  per 
bit,  normalized  to  the  noise  power,  for  different  numbers  of  antennas  per  terminal.  For  a  fair  comparison, 
all  curves  refer  to  the  same  overall  transmitted  power.  The  constellation  is  QPSK.  The  overall  number 
of  transmit/receive  antennas  is  also  the  same  in  all  cases.  More  specifically,  we  used  the  following 
combinations:  n  =  100  and  ny  =  1  (dotted  line),  n  =  50  and  nr  =  2  (dashed  line),  and  n  =  25  and 


^D^{x)  denotes  the  inverse  of  Dxix)  and  it  certainly  exists  because  in  this  case  Dx{x)  is  strictly  monotone. 
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Fig.  2.  Node  density  (number  of  nodes  per  area  unit)  versus  Si, /No,  for  different  number  of  antennas  per  terminal. 


Fig.  3.  Node  density  (number  of  nodes  per  area  unit)  versus  Si /No,  for  different  number  of  antennas  per  terminal  (nr  =  1 
or  2)  and  different  transmission  rates  (  M  =  4, 16, 64). 


riT  =  4  (solid  line).  The  conneetivity  is  insured  with  prohahility  pcon  =  0.99  and  the  out-of-service 
event  refers  to  a  target  BER  of  10“^  and  it  is  required  to  occur  with  a  maximum  time  percentage  of 
Pout  =  10“^.  We  can  see,  from  Figure  2,  the  advantage  of  diversity  that  makes  possible  a  considerable 
decrease  of  the  nodes  density.  Clearly,  terminals  with  multiple  antennas  provide  a  larger  coverage  than 
single  antenna  terminals,  but  at  the  cost  of  increased  complexity. 

The  connectivity  is  also  a  function  of  the  bit  rate.  As  an  example,  in  Figure  3  we  show  the  minimum 
SNR  required  for  connectivity,  assuming  an  efficiency  of  2,4,  and  6  bits/sec/Hz,  achieved  using  4,  16, 
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or  64-QAM  constellations.  As  expected,  an  increase  of  rate  requires  an  increase  of  node  density,  for  a 
given  power  budget.  Hence,  the  bit  rate  may  result  as  a  compromise  between  the  number  of  antennas, 
node  density,  energy  per  node,  and  complexity. 


E.  Cooperative  Communications 

The  next  question  is  “What  can  we  do  to  improve  the  connectivity  if  we  have  only  single  antenna 
transceivers?”.  We  can  resort  to  cooperation  among  nearby  terminals.  The  idea  is  pictorially  sketched  in 
Figure  4,  where  we  see  three  nodes.  A,  B,  and  C.  In  the  absence  of  any  cooperation,  each  node  covers 
a  circle  of  radius  ro  given  by  (13).  For  example,  in  Figure  4,  A  and  B  are  connected,  but  C  is  isolated. 
However,  if  nodes  A  and  B  cooperate,  they  can  give  rise  to  a  virtual  transmit  array  capable  of  covering 


a)  coverage  without  cooperation 


Fig.  4.  Coverage  in  cooperating  networks. 

an  area  larger  than  that  covered  with  a  single  antenna  [22]  (see  also  [50]  and  the  references  therein).  The 
idea  is  represented,  pictorially,  in  Figure  4  (b),  where  the  bigger  circle  is  the  area  covered  by  a  system 
located  in  the  center  of  gravity  of  the  nodes  A  and  B,  with  a  bigger  radius  resulting  from  the  use  of  a 
MISO  system  with  two  transmit  and  one  receive  antenna.  Proceeding  as  in  the  previous  section,  denoting 
with  nreiay  the  number  of  cooperating  (relay)  nodes,  we  have  the  potential  of  diversity  gain  nreiay  + 1  (the 
relays  plus  the  source  itself).  The  existence  of  the  bigger  circle  is  a  result  of  the  cooperation  between  A 
and  B.  Thanks  to  cooperation,  a  disconnected  network  may  become  connected,  as  shown  in  the  example 
of  Figure  4  (b),  using  the  same  overall  transmit  power. 

Let  us  now  quantify  how  much  is  the  coverage  increase,  due  to  cooperation,  and  how  this  affects  the 
connectivity.  If  we  have  rireiay  relays  cooperating  with  a  source  and  a  destination  with  one  receive 
antenna,  the  (maximum)  diversity  gain  is  ut  '■=  n^eXay  +  1-  Repeating  derivations  similar  to  the  ones 
described  in  the  previous  section,  with  the  only  exception  that  now  we  only  have  transmit  diversity,  but 
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no  receive  diversity  because  each  receiver  has  a  single  antenna,  the  coverage  radius  in  case  of  cooperation 
is 

=  f3r  cov,  (25) 


gnSb 


Nout  log{cM/Po) 


{ut'-Pc 


out ) 


where  Vcov  is  given  by  (14)  and 


(3:= 


(nr!) 


p(nT-l)/nT 


(26) 


Therefore,  the  coverage  increases  by  a  factor  (3  that  depends  on  the  number  of  cooperating  nodes  and 


on  the  desired  out-of-service  probability. 


The  effect  of  the  coverage  increase  on  the  network  connectivity  is  illustrated  in  Figure  5,  where  we  report 
the  connection  probabilities  obtained  without  cooperation  (dashed  line)  and  with  cooperation  (solid  line). 
In  case  of  cooperation,  we  considered  only  the  case  of  no  more  than  two  cooperating  terminals.  The 
probabilities  shown  in  Figure  5  have  been  estimated  over  a  set  of  200  independent  network  realizations. 
As  a  comparison  term,  we  report  in  Figure  5  the  connection  probability  of  a  non-cooperative  network, 
but  having  a  coverage  radius  /3rcov  (dotted  line).  We  can  see  that  cooperation  between  pairs  of  radio 


Fig.  5.  Connection  probabilities  versus  SNR  (dB):  without  cooperation  (dashed  line),  with  cooperation  (solid  line),  and  without 
cooperation,  but  using  /3rcov  instead  of  rcov  (dotted  line). 


nodes  is  sufficient  to  yield  an  SNR  gain  of  approximately  seven  dB.  Interestingly,  the  curve  obtained 
without  cooperation,  but  with  a  coverage  radius  /3rcov  has  approximately  the  same  connectivity  as  the 
cooperative  case,  where  the  coverage  of  each  node  is  rcov-  Recalling,  from  (25),  that  the  transmitted 
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power  is  proportional  to  r^,  if  r  is  the  coverage,  we  infer  that  cooperation  among  pairs  of  terminals 
yields  an  improvement  in  terms  of  transmitted  power  approximately  equal  to  /3^,  when  nr  =  2. 

In  conclusion,  the  approximate  expressions  derived  in  this  paper  have  been  shown  to  fit  quite  well 
the  behavior  of  wireless  networks  having  a  random  topology,  modeled  as  a  2D  Poisson  point  process, 
in  terms  of  connectivity.  The  formulas,  albeit  approximate,  help  to  predict  the  gain  achievable  by  using 
multiantenna  terminals,  either  real  or  virtual,  i.e.  through  cooperative  communications.  We  have  derived 
expressions  for  the  density  with  which  the  network  nodes  should  be  deployed  in  order  to  guarantee  the 
network  connectivity  with  a  desired  probability.  In  this  paper,  we  have  only  shown  results  concerning 
Rayleigh  fading  channels.  We  have  extended  the  analysis  to  Nakagami-m-fading  channels  and  we  showed 
that  the  advantage  decreases  as  the  index  m  increases,  i.e.  as  the  channel  tends  to  be  less  and  less  random. 


F.  Cooperation  gain 


The  probability  that  k  nodes  fall  within  any  given  circle  of  radius  tq  is  then 


(27) 


For  any  given  transmit  power,  (27)  gives  the  probability  of  finding  k  relays.  The  coverage  radius  rg,  on 
ifs  furn,  depends  on  fhe  energy  Sb,  as  well  as  on  fhe  fransmission  rafe  (fhroughouf  fhe  fwo  coefficienfs 
Cm  and  ^m)-  Clearly,  increasing  8b,  the  coverage  increases,  but  more  energy  is  wasted  only  to  send 
information  towards  the  relays  instead  of  the  final  infended  desfinafion.  There  is  in  general  an  opfimal 
energy  disfribufion  befween  fhe  fwo  phases  where  fhe  source  sends  dafa  fo  fhe  relay  and  when  source 
and  relays  fransmif  fogefher  fowards  fhe  desfinafion.  Lef  us  denofe  by  8t  the  total  energy  per  bit  for  all 
cooperating  nodes.  Introducing  the  real  coefficient  (3,  with  0  <  /?  <  1,  we  indicate  with  PSt  the  portion 
of  the  total  energy  dedicated  to  send  information  from  the  set  of  cooperating  nodes  to  the  destination 
and  with  (1  —  P)St  the  energy  spent  by  the  source  to  send  data  to  the  relays. 


If  k  is  the  number  of  relays  that  receive  the  data  from  the  source  with  the  prescribed  reliability,  the  k 
relays  plus  the  source  can  then  transmit  together  towards  the  destination  as  if  they  were  the  transmit 
antennas  of  a  single  user^.  Denoting  with  hi  the  channel  coefficients  from  source  and  relays  towards  the 
destination,  using  a  full-diversity  space-time  coding  scheme,  such  as,  e.g.,  orthogonal  coding,  the  error 

^  Since  a  node  is  chosen  as  a  relay  only  if  its  BER  is  below  a  given  threshold,  with  a  given  outage  probability,  we  assume 
here  that  the  errors  at  the  relay  nodes  are  negligible. 
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probability  at  the  receiver  is^ 

Pe{k  +  l]h)  =  cQ 


P£t 


(28) 


Assuming  that  the  channels  are  statistically  independent,  the  expected  value  of  Pf.{k  +  1)  :=  Et^{Pe{k  + 
1;  h)}  is  [33] 

aVM-1  /I ^  (  k  + 


Pe{k  +  1)  = 


\/Mlog2(M) 


■ 


E 

m=0  \  ^ 


m 


1  +  /i 


where 


9  ■= 


3P£Tlog2{M)al 


2>(3£T\og2{M)(Tl  +  2{M  -  l){k  + 


The  final  average  error  probability  at  the  destination,  in  case  of  cooperation,  is  then 

CXD 

Pe  =  '^Pro{k)Pe{kPl), 


(29) 


(30) 


(31) 


k=0 


where  Pro{k)  is  given  by  (27),  with  vq  given  by  (13)  setting  Sb  =  {1  —  P)£t,  whereas  Pe{k  + 1)  is  given 
by  (29).  An  example  of  average  BER  is  reported  in  Fig.  6,  for  a  network  of  nodes  with  different  node 


Fig.  6.  Average  BER  at  the  final  destination  in  a  cooperating  network  for  different  values  of  node  density. 


''in  case  of  coordinated  transmission  from  fc  +  1  nodes,  we  normalize  the  transmit  power  of  each  node  by  fe  +  1,  so  that  the 


overall  radiated  power  is  independent  of  k. 
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densities.  We  can  clearly  see  the  advantage  of  using  cooperation  with  respect  to  the  non-cooperative  case 
and  how  the  gain  increases  as  the  node  density  increases.  At  high  SNR,  the  dominant  term  in  (31)  is  the 
term  with  /c  =  0,  as  it  corresponds  to  the  slowest  decay  rate  of  the  average  hit  error  rate.  The  term  Pe{^) 
goes  like  1/ SNR,  at  high  SNR.  Hence,  in  a  random  network  there  is  no  real  diversity  gain.  Nevertheless, 
since  Pe(l)  is  multiplied  hy  Pi?(0),  there  still  is  a  coding  gain  equal  to  1/pr{^),  i.e. 

Gc  =  e^P^o.  (32) 

Hence,  cooperation  introduces  diversity,  in  the  sense  of  the  existence  of  multiple  paths  for  the  transmitted 
data,  hut  unlike  conventional  systems,  this  results  in  a  coding  gain,  not  in  a  diversity  gain.  The  coding 
gain  is  always  greater  than  one  and  it  grows  exponentially  with  the  increase  of  the  relay  nodes  density. 
To  obtain  a  higher  Gc,  for  a  given  p,  it  is  necessary  to  increase  the  coverage  radius  tq.  This  requires  that 
more  energy  is  used  in  the  first  phase,  when  S  sends  data  to  the  relays.  If  there  is  a  constraint  on  the  total 
energy  St,  it  is  interesting  to  see  what  is  the  optimal  power  distrihution  between  the  two  transmission 
phases,  acting  on  /3. 

As  an  example,  in  Fig.  7  we  show  the  average  BER,  for  a  given  SNR  at  the  destination,  as  a  function 
of  j3.  We  can  see  that,  depending  on  the  final  SNR,  there  is  an  optimal  /3. 


Fig.  7.  Average  BER  as  a  function  of  /3,  for  different  SNR  values. 
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II.  Distributed  space-time  coding 

The  choice  of  the  most  appropriate  space-time  coding  technique  to  he  used  in  a  relay  network  depends 
on  several  factors.  As  with  conventional  space-time  coding,  the  choice  depends  on  the  desired  trade¬ 
off  between  rate,  diversity  and  receiver  complexity.  The  number  nR  of  real  receive  antennas  plays 
a  fundamental  role.  If  ur  =  1,  the  cooperation  induces  a  virtual  multi-input/single  output  (MISO) 
communication,  introducing  diversity  into  the  system,  but  not  a  rate  gain.  Conversely,  if  >  1,  we  can 
think  of  schemes  capable  of  both  diversity  and  rate  gains,  exploiting  the  resulting  virtual  MIMO  structure. 
Basically,  one  could  choose  among  the  following  classes  of  STC  techniques:  i)  Orthogonal  STC  (OSTC) 
[7],  as  a  strategy  that  maximizes  the  diversity  gain  and  it  minimizes  the  receiver  complexity;  ii)  full- 
rate/full  diversity^  codes  (FRFD)  [34],  [35],  as  codes  that  yield  maximum  diversity  gain  and  transmission 
rate,  but  with  high  receiver  complexity;  iii)  V-BLAST  codes  [37],  as  a  technique  that  maximizes  the 
rate,  sacrificing  the  diversity  gain,  but  with  limited  receiver  complexity.  Alternatively,  one  could  use  the 
trace-orthogonal  design  [38]  as  a  flexible  way  to  trade  complexity,  bit  rate  and  bit  error  rate.  It  is  worth 
noticing  that  the  optimal  trade-off  among  these  alternative  strategies,  in  the  distributed  case  does  not 
coincide,  necessarily,  with  the  trade-off  achievable  with  conventional  space-time  coding. 

The  first  evident  difference  between  conventional  STC  and  DSTC  is  the  presence  of  the  time  slot  necessary 
for  the  exchange  of  data  between  source  and  relays.  This  induces  an  inevitable  rate  loss.  To  reduce  this 
loss,  it  is  necessary  to  allow  for  the  re-use  of  the  same  time  slot  by  more  than  one  set  of  source -relay 
pairs.  In  Fig.  8,  we  show,  as  an  example,  three  sources  (circles)  and  some  potential  relays  (dots).  If  the 
relays  are  associated  to  the  nearest  sources®  and  the  sources  are  sufficiently  far  apart,  we  can  assign  the 
same  time  slot  for  the  exchange  of  information  between  each  source  and  its  own  relay.  Clearly,  this  does 
not  prevent  the  interference  between  different  source/relay  pairs.  In  general,  the  relay  discovery  phase 
should  follow  a  strategy  that  gives  rise  to  many  spatially  separated  micro-cells,  as  in  Fig.  8,  where  each 
source  acts  as  a  local  base  station  broadcasting  to  its  relays,  who  may  get  interference  from  other  micro¬ 
cells  (sources).  The  need  to  limit  the  coverage  of  each  source,  in  its  relay  discovery  phase,  is  also  useful 
because:  i)  less  power  is  wasted  in  the  source-relay  link;  ii)  there  are  less  synchronization  problems  in 

^The  term  full  rate  here  is  used  in  the  same  sense  as  [34],  [35]  and  it  means  that  a  transmitter  with  tit  antennas  transmits 
rty  symbols  in  ut  time  slots.  This  does  not  imply  anything  about  the  final  BER  and  it  has  then  to  be  distinguished  by  the 
information  rate  concept  used  in  [36]. 

®We  will  comment  later  on  the  meaning  of  distance  between  source  and  relay,  as  it  has  to  take  into  account  also  the  channel 
fading. 
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Fig.  8.  Geometry  of  cooperative  network  with  three  source  nodes  (circles),  potential  relays  (dots)  and  one  destination  (D). 


the  final  link  towards  the  destination;  iii)  there  is  less  interference  between  different  source/relay  sets. 


We  now  compute  the  rate  loss  for  different  cooperation  structures,  assuming  that  there  are  N  simultaneous 
source/relays  pairs  sharing  the  same  time  slot.  We  denote  with  rs2R  and  Tsr2D  the  duration  of  the  S2R 
and  SR2D  time  slots.  Tg  is  the  symbol  duration  in  all  slots.  For  a  given  bit  rate,  the  durations  depend 
on  the  constellation  order  used  in  the  different  slots.  We  denote  with  Q  and  M  the  constellation  orders 
used  in  the  S2R  and  in  the  SR2D  slots,  respectively.  We  consider  a  frame  containing  both  S2R  and  the 
N  SR2D  links  has  then  a  duration  Tp  =  Ts2R  +  NTsr2D-  The  rate  reduction  factor,  with  respect  to  the 
non-cooperative  case,  in  a  TDMA  context,  is  then 

-/VrsR2D 


7?  = 


(33) 


Ts2R  +  -^Tsr2D 

Clearly,  the  rate  loss  can  be  reduced  by  decreasing  Ts2R,  i-e.  by  increasing  Q,  or  by  increasing  N.  In  the 
first  case,  the  relay  needs  a  higher  Signal-to-Noise-plus-Interference-Ratio  (SNIR).  In  the  second  case, 
SNIR  decreases  at  the  relay,  as  there  is  more  interference.  In  both  cases,  it  is  less  likely  to  discover 
a  relay  with  sufficient  SNIR.  Hence,  the  right  choice  has  to  result  from  a  trade-off  between  rate  and 
performance.  We  discuss  now  in  detail  the  alternative  DSTC  strategies,  corresponding  to  OSTC,  FRFD 
and  V-BLAST  codes.  In  all  cases,  we  denote  with  s{n)  the  sequence  of  symbols  sent  by  S  during  the 
SR2D  slot,  whereas  s{n)  indicates  the  estimate  of  s{n)  performed  at  the  relay.  For  simplicity,  we  refer 
to  a  context  where  each  source  cooperates  with  only  one  relay. 
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A.  Distributed  Orthogonal  STC  (D-OSTC) 

D-OSTC  guarantees  maximum  receiver  simplicity  and  full  diversity  and  it  can  be  implemented  also 
when  the  final  destination  has  a  single  antenna.  D-OSTC  was  proposed  in  [22],  [19],  [25],  where  the 
relays  are  essentially  error-free,  and  in  [29],  where  the  decoding  errors  at  the  relay  are  explicitly  taken 
into  account.  D-OSTC  transmits  2  symbols  over  two  successive  time  periods,  so  that  Tsr2d  =  2Ts  and 
Ts2R  =  2  log2(M)rs/log2((5).  The  sequence  transmitted  by  the  source-relay  pair  is 

s(n)  -s*(n+l) 

.  (34) 

s{n  -h  1)  s*{n) 

The  first  row  of  this  matrix  contains  the  symbols  transmitted  by  the  source,  whereas  the  second  row 
refers  to  the  symbols  transmitted  by  the  relay  (different  columns  refer  to  successive  time  instants).  The 
overall  bit  rate,  incorporating  also  the  rate  loss,  is 

R  =  - Ti  b/s/Hz.  (35) 

2N  -|-  2  log2  M/  log2  Q 

B.  Distributed  full  rate/full  diversity  (D-FRFD) 

If  the  final  destination  has  2  antennas,  there  is  a  virtual  2x2  MIMO,  with  the  possibility  of  increasing 
the  rate.  This  can  be  achieved,  for  example,  using  distributed-FRFD  or  distributed-BLAST.  With  D-FRFD, 
the  pair  S-R  transmits  4  symbols  over  two  consecutive  time  periods.  The  transmitted  matrix  is  [34]  or 


s(n)  +  (/9s(n  +  1)  6*  (s(n  +  2)  +  (/3s(n  +  3)) 

,  (36) 

0  (s(n -h  2)  —  (/9s(n  +  3))  s{n)  —  ifs{n  +  l) 

where  ip  =  9  =  are  two  rotation  parameters  (see,  e.g.  [34]  or  [35],  for  the  choice  of  p  and  9). 


The  bit  rate  is 


D  -  4A^log2M 

2iV  +  41og2M/log2Q 


C.  Distributed  BLAST  (D-BLAST) 


We  consider  here  the  version  of  BLAST  where  two  independent  streams  of  data  are  transmitted  from 
the  two  antennas.  In  its  distributed  version,  D-BLAST  requires  that  the  relay  receives  only  half  of  the 
bits  to  be  transmitted.  This  implies  an  advantage  with  respect  to  D-FRFD,  as  it  allows  us  to  reduce 
the  duration  of  the  S2R  time  slot.  The  price  paid  with  respect  to  D-FRFD  is  that  D-BLAST  is  not  full 


diversity.  The  transmitted  matrix  in  the  D-BLAST  case,  is 

s{n)  s{n  +  2) 
s(n  +  l)  s{n  +  3) 


(38) 
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and  the  bit  rate  is 


R  = 


4N  log2  M 


b/s/Hz. 


(39) 


2A^  +  21og2  M/  log2Q 

Comparing  the  transmission  rates  of  all  the  distributed  schemes,  for  a  given  choice  of  the  constellation 
orders  Q  and  M,  we  see  that  D-BLAST  has  the  highest  transmission  rate. 


D.  Relay  discovery  phase 

This  section  describes  a  possible  resource  discovery  phase.  A  source  looking  for  potential  relays  starts 
sending  a  sounding  signal  to  verify  whether  there  are  available  neighbors.  The  sounding  signal  is  a 
pseudo-noise  code  identifying  the  source.  A  potential  relay  may  receive  the  sounding  signals  from  more 
than  one  source.  The  radio  nodes  available  to  act  as  relays  compute  the  signal-to-noise  plus  interference 
ratio  (SNIR)  for  each  source  This  step  requires  the  node  to  be  able  to  separate  the  signals  coming  from 
different  sources.  This  is  made  possible  by  the  use  of  orthogonal  codes.  The  potential  relays  retransmit  an 
acknowledgment  signal  back  only  to  those  sources  whose  SNIR  exceeds  a  certain  threshold.  The  source 
receives  then  the  acknowledgments  and  the  relative  SNIR  from  all  potential  relays  and  it  decides  which 
relays  to  use.  This  phase  insures  that  the  relay,  once  chosen,  is  sufficiently  reliable.  Given  the  variability 
of  the  wireless  channel,  this  operation  has  to  be  repeated  at  least  once  every  channel  coherence  time. 
To  avoid  excessive  complications,  a  node  may  act  as  a  relay  for  no  more  than  one  source.  The  basic 
philosophy  we  follow  to  discover  relays  is  that  source  and  relays  should  be  as  close  as  possible.  This 
is  justified  by  the  following  concurring  reasons:  i)  less  power  is  wasted  in  the  S2R  slot;  ii)  there  are 
less  synchronization  problems  in  the  final  SR2D  slot;  iii)  there  is  less  interference  between  different 
source/relay  sets.  In  summary,  this  relay  discovery  phase  creates  many  spatially  separated  micro-cells 
where  each  source  acts  as  a  local  base  station  broadcasting  to  its  relays,  who  may  get  interference  from 
other  micro-cells  (sources). 

^The  SNIR  may  be  evaluated  for  each  channel  realization  or  in  average  sense,  considering  the  channel  statistics  such  as  mean 
and  covariance.  The  first  option  provides  better  performance,  but  it  requires  more  frequent  channel  estimation  updates  than  the 
second  option. 
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III.  Regenerative  relaying  with  block  distributed  orthogonal  space-time  coding  in 

SINGLE-USER  SCENARIO 

In  this  section,  we  analyze  in  detail  an  example  of  Distributed  Space-time  Coding,  using  D-OSTC  and 
regenerative  relay  [29].  We  assume  that  the  source  has  found  its  own  relay^  and  that  the  interference  due 
to  the  other  active  source-relay  links  is  negligible.  We  consider,  thus,  a  two-hop  relay  channel,  composed 
of  a  source  (S),  a  relay  (R)  and  a  destination  (D). 

Differently  from  common  STC,  with  DSTC,  i)  regenerative  relays  might  make  decision  errors,  so  that 
the  symbols  transmitted  from  R  could  be  affected  by  errors;  ii)  the  links  between  S  and  D  and  between 
R  and  D  do  not  have  the  same  statistical  properties,  in  general;  iii)  even  if  S  and  R  are  synchronous, 
their  packets  might  arrive  at  D  at  different  times,  as  S  and  R  are  not  co-located.  In  the  following,  we 
will  address  all  these  problems  specifically. 

We  illustrate  the  proposed  transmission  protocol  by  referring  to  a  TDD  scheme,  but  the  same  consid¬ 
erations  could  apply  to  an  FDD  mode.  In  a  TDD  system,  each  frame  is  subdivided  in  consecutive  time 
slots:  In  the  first  slot  S  transmits  and  R  receives;  in  the  second  slot,  S  and  R  transmit  simultaneously. 
We  describe  the  distributed  space-time  protocol  within  the  following  setup:  (al)  all  channels  are  FIR 
of  (maximum)  order  and  time-invariant  over  at  least  a  pair  of  consecutive  blocks;  (a2)  the  channel 
coefficients  are  i.i.d.  complex  Gaussian  random  variables  with  zero  mean  and  variance  1/d?,  where  d 
is  the  link  length;  (a3)  the  information  symbols  are  i.i.d.  BPSK  symbols  that  may  assume  the  values 
A  OT  —A  with  equal  probability^;  (a4)  the  received  data  are  degraded  by  additive  white  Gaussian  noise 
(AWGN);  (a5)  the  channels  are  perfectly  known  at  the  receive  side  and  are  unknown  at  the  transmit 
side;  (a6)  the  transmission  scheme  for  all  terminals  is  blockwise,  where  each  block  is  composed  of  M 
symbols,  incorporating  a  cyclic  prefix  of  length  L  equal  to  the  sum  of  the  relative  delay  with  which 
packets  from  S  and  R  arrive  at  D  plus  the  maximum  channel  order  Lh- 

We  will  use  the  following  notation.  We  denote  with  hgd,  hgr,  and  hrd,  the  impulse  responses  between  S 
and  D,  S  and  R,  and  R  and  D,  respectively.  Each  block  of  symbols  s{i)  has  size  M  and  it  is  linearly 
encoded,  so  as  to  generate  the  A^-size  vector  Xs(n)  :=  Fs(n),  where  F  is  the  N  x  M  precoding  matrix. 
A  CP  of  length  L  >  L/j  is  inserted  at  the  beginning  of  each  block,  to  facilitate  elimination  of  inter-block 
interference,  synchronization,  and  channel  equalization  at  the  receiver.  denotes  the  pseudo-inverse  of 

^Thanks  to  the  assumption  that  the  relay  nodes  are  not  error-free,  the  probability  of  finding  a  relay  in  the  discovery  phase 
increases,  with  respect  to  the  schemes  that  use  a  relay  only  if  it  is  error-free. 

^Assumption  (a3)  is  made  only  for  simplifying  our  derivations,  but  there  is  no  restriction  to  use  higher  order  constellations. 
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A;  indicates  the  real  part  of  x;  when  applied  to  a  vector,  is  the  vector  whose  entries  are  the 

real  part  of  the  entries  of  x. 

D-OSTC,  for  frequency  selective  channels,  works  as  follows.  During  the  first  time  slot,  S  sends, 
consecutively,  the  two  A^-size  information  symbols  blocks  s{i)  and  s{i  +  1).  The  blocks  are  linearly 
encoded  using  the  precoding  matrix  F,  so  that  the  corresponding  transmitted  blocks  are  Xs{n)  :=  Fss{n), 
with  n  =  i,i  +  1.  Under  (a6),  after  removing  the  guard  interval  at  the  receiver,  the  A^-size  vectors  ^^(n) 
received  from  R  are 

yr{n)  =  HsrFss{n)  +  Wr{n),  n  =  i,i  +  l,  (40) 

where  Wr{n)  is  the  additive  noise  at  the  relay.  Thanks  to  the  insertion  of  the  CP,  the  channel  matrix 
Hsr  is  A^  X  A^  circulant  Toeplitz  and  it  is  diagonalized  as  Hgr  =  W AgrW^ ,  where  W  is  the  N  x  N 
IFFT  matrix  with  {W}ki  =  whereas  is  the  N  X  N  diagonal  matrix,  whose  entries 

are  =  Ef=o“' ^-(0 

The  relay  node  decodes  the  received  vectors  and  provides  the  estimated  vectors  s{i)  and  s{i  +  1). 

During  the  successive  time-slot,  S  and  R  transmit  simultaneously,  using  a  block  Alamouti’s  strategy 
[39].  More  specifically,  in  fhe  firsl  half  of  fhe  second  lime  slol,  S  Iransmits  Xs{i  +  2)  =  ai  Fs{i)  and 
R  Iransmits  Xr{i  +  2)  =  a2  Fs{i  +  1).  In  the  second  half,  S  transmits  Xs{i  +  3)  =  ai  Gs*{i  +  1)  while 
R  transmits  Xr{i  +  3)  =  —a2  Gs*{i).  To  guarantee  maximum  spatial  diversity,  the  two  matrices  G  and 
F  are  related  to  each  other  by  G  =  Ji^*,  as  in  [39],  where  J  is  a  time  reversal  (plus  a  one  chip  cyclic 
shift)  matrix.  If  N  is  even  J  has  all  null  entries  except  the  elements  of  position  (1, 1)  and  {k,N  —  k  +  2), 
with  k  =  2, ...  ,N,  which  are  equal  to  one.  If  N  is  odd,  J  is  the  anti-diagonal  matrix.  The  two  real 

coefficients  ai  and  a2  are  related  to  each  other  hy  a\  +  a\  =  1.  They  are  introduced  in  order  to  have 

a  degree  of  freedom  in  the  power  distribution  between  S  and  R,  under  a  given  total  transmit  power.  In 
Section  III-C,  we  will  show  how  to  choose  ai  (and  then  02)  in  order  to  minimize  the  final  average  BER. 

Affer  discarding  fhe  CP  and  using  (al),  the  blocks  received  by  D  in  the  two  consecutive  time-slots 
f  -|-  2,  and  f  -|-  3  are  given  by 

Vdi'^  +  2)  =  aiH sdF s{i)  +  a2H rdF s{i  -|-  1)  -|-  Wd{i  +  2) 
yd(*  +  3)  =  aiHsdGs*{i  + 1)  -  a2HrdGs*{i)  +  Wd{i +  3), 
where  Hgd  and  Hj.d  refer  to  the  channels  between  S  and  D  and  between  R  and  D,  respectively. 
Exploiting,  again,  the  diagonalizations  Hgd  =  W AgdW^  and  Hrd  =  W ArdW^ ,  if  we  pre-multiply 
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in  (41)  y^{i  +  2)  by  and  y*^{i  +  3)  by  W'^ ,  we  get 


W^y^{i  +  2)  =  aiAsdFs{i)  +  a2ArdFs{i  +  1)  +  W^Wd{i  +  2) 

W^ylii  +  3)  =  aiAlaG*s{i  +  1)  -  a2A;aG*s{i)  +  +  3), 

where  F  :=  F  and  G  :=  W^G.  For  the  sake  of  simplicity,  we  assume  that  OFDM  is  performed 

at  both  S  and  R  nodes,  so  that  N  =  M,  F  =  W  and  thus  F  =  I]\f  and  G  =  W.  We  also  introduce 

the  orthogonal  matrix 

OtlAgd  Ot2Ayd 

such  that  A^A  :=  120 A^,  where  A^  :=  a\\Asd\‘^ +oi\\Ard\^ ,  whereas  0  denotes  the  Kronecker  product. 
We  introduce  also  the  unitary  matrix*^  Q  :=  A  {I2  0  A  ^),  satisfying  the  relationships  Q^Q  =  I2N 
and  A  =  I2  0  A.  Exploiting  the  above  equalities  and  multiplying  the  vector  u  :=  [{W^y^{i  + 
2))^,  {W'^y*^{i  +  3))^]^  by  the  matrix  ,  without  compromising  the  decision  optimality  (because  of 
the  unitarity  of  Q),  we  get 


r{i  +  1) 


~  jf:  ~ 

~^sd^rd 

\~Ard? 

~  jf:  ~ 

■^sd-^rd 

s  + 

1 

~AsdA^^ 

- 1 

(43) 


where  s  :=  [s(i)'^,s(i  +  1)'^]'^,  s  :=  [s(i)^,s(i  +  1)'^]^,  A^d  ■=  aiA^^A  A^d  :=  a2A^rfA 
w  :=  [w'^{i),  w^{i  +  1)]"^  =  +  2),  w^{i  +  3)]'^.  As  expected,  the  previous  equations  reduce 

to  the  classical  block  Alamouti  equations,  see  e.g.  [39],  the  two  transmit  antennas  use  the  same  power, 
i.e.,  ai  =  a2,  and  there  are  no  decision  errors  at  the  relay  node,  i.e.,  s(n)  =  s(n),n  =  i,i  +  1. 

Since  is  unitary,  if  w  is  white,  w  is  also  white,  with  covariance  matrix  C^,  =  a^l2N-  Furthermore, 
since  all  matrices  A  appearing  in  (43)  are  diagonal,  the  system  (43)  of  2N  equations  can  be  decoupled 
into  N  independent  systems  of  two  equations  in  two  unknowns,  each  equation  referring  to  a  single  sub¬ 
carrier.  More  specifically,  introducing  the  vectors  :=  [rk{i) ,  rk{i+l)]'^ ,  Sk  :=  [sk{i),  Sk{i+l)]'^,  Sk  := 
[sfc(i),  Sfc(i-|- 1)]'^  and  Wk  :=  [wk{i),  Wk{i  +  1)]'^,  referring  to  the  fc-th  sub-carrier,  with  k  =  0, . . . ,  N  —  1 
(for  simplicity  of  notation,  we  drop  the  block  index  and  we  set  A^^  =  Asd{k,  k)  and  A^d  =  ^rd{k,  k)). 


'°We  suppose  that  the  channels  do  not  share  common  zeros  on  the  grid  Zq  =  q  integer,  so  that  A  is  invertible. 


December  5,  2006 


DRAFT 


24 


(43)  is  equivalent  to  the  following  systems  of  equations 

Since  the  noise  vector  Wk  is  also  white  with  covariance  matrix  =  (t‘^I2N,  and  there  is  no  inter- 
symhol  interference  (ISI)  between  vectors  and  corresponding  to  different  suh-carriers,  rj.  represents 
a  sufficient  statistic  for  the  decision  on  the  transmitted  symbols  vector 


+ 


A!  A 


sd^^rd 


|A, 


H  +  Wk 


(44) 


rd\ 


/  |A,rf|2  -klJ^rd  \ 
\  KdKd  lA.dP  ) 


A.  ML  detector 


We  derive  now  the  structure  of  the  maximum  likelihood  (ML)  detector  at  the  final  destinafion.  Besides 
fhe  previous  assumpfions,  we  assume  also  fhat  D  has  perfecf  knowledge  of  fhe  vector  of  error  probabilifies 
Pei{k)  and  pe^{k),  k  =  .  ,N  —  1,  occurring  af  fhe  relay.  This  requires  an  exchange  of  informalion 

befween  R  and  D.  This  informafion  has  fo  be  updafed  wifh  a  rale  depending  on  fhe  channel  coherence 
lime.  Later  on,  we  will  show  an  alternative  (sub-optimum)  detection  scheme  that  does  not  require  such 
a  knowledge. 

We  denote  with  A  the  set  of  all  possible  transmitted  vectors  and  with  Pei{k)  and  Pe2{k)  the 
conditional  (to  a  given  channel  realization)  error  probabilities,  at  the  relay  node,  on  Sfc(l)  and  5^(2), 
respectively.  After  detection,  at  the  node  R,  we  have  Sk{l)  =  Sfc((),  with  probability  (1  —  pei{k)),  or 
Sk{l)  =  —Sk{l),  with  probability  Pei{k),  I  =  1,2.  Since  the  symbols  are  independent,  the  probability 
density  function  of  the  received  vector  z,  conditioned  to  having  transmitted  s^,  is  [40] 

fz\S,{z\Sk)=  [{1  -  Pel{k)){l  -  Pe2{k))  exp  {  - 1 Z  -  Afc(l ,  1)  Sfc  |  } 

TT 

+Pel{k)pe2{k)  exp  {-|z  -  Afc(-1, -l)Sk\^/ol} 

+(1  -pei(A:))pe2(fe)exp{-|z  -  Afc(l,  -l)sk\^/al} 

+Pei{k){l  -  Pe2{k))  exp  {-\z-  Afc(-1,  l)sk\^/al]]  , 

(45) 


where  Ak{9i,92)  is  defined  as  follows 


Aa:(0i,02)  = 


|AsrfP  +  \Ard\^9i,  A*J^rd92  “  ^Id^rd 
AsdKd  -  AsdKd^i,  |A,dp  +  \Kd?92 
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Based  on  (45),  the  ML  detector  is 

Sfc  =  arg  max  { /^^|s^(rfc|sfc)}  .  (46) 

Note  that,  thanks  to  the  orthogonal  space-time  block  coding  strategy,  the  optimal  detector  preserves  the 
receiver’s  simplicity,  because,  under  (al)-(a6),  the  ML  solution  performs  an  exhaustive  search  only  among 
four  possible  transmitted  vectors  s^’s. 

B.  Sub-optimum  detector 

The  ML  detector  described  above  requires  the  knowledge,  at  the  destination  node,  of  the  set  of  error 
probabilities  Pei{k)  and  Pe2{k),  with  A:  =  0, . . . ,  —  1.  If  this  knowledge  is  not  available,  a  sub-optimum 

scalar  detector  can  be  implemented,  instead  of  the  ML  detector.  More  specifically,  the  decision  on  the 
transmitted  symbol  Sk{n)  can  be  simply  obtained  as 

s(n)  =  sign  {3f?[r(n)]}  ,  n  =  i,i-\-l,  (47) 

where  r(n)  is  given  by  (43).  Note  that,  for  high  SNR  at  the  relay  (i.e.  when  R  makes  no  decision 
errors),  the  symbol-by-symbol  decision  in  D  becomes  optimal  and,  thus,  the  decoding  rule  (47)  provides 
the  same  performance  as  the  optimal  receiver  (46).  When  the  decision  errors  at  the  relay  side  cannot 
be  neglected,  the  sub-optimal  receiver  introduces  a  floor  in  the  bit-error-rate  (BER)  curve,  because  the 
symbol-by-symbol  decision  (47)  treats  the  wrong  received  symbols  as  interference.  The  choice  between 
the  decoding  rules  (46)  and  (47)  should  then  result  as  a  trade-off  between  performance  and  computational 
complexity,  taking  into  account  the  need,  for  the  ML  detector,  to  make  available,  at  the  destination  node, 
the  error  probabilities  of  the  relay  node.  We  will  show  a  comparison  between  ML  and  sub-optimum 
strategies  in  Section  III-D. 

C.  Power  Allocation  between  Source  and  Relays 

While  in  conventional  STC,  the  transmit  antennas  typically  use  the  same  power  over  all  the  transmit 
antennas,  with  DSTC  it  is  useful  to  distribute  the  available  power  between  source  and  relay  as  a  function 
of  their  relative  position  with  respect  to  the  final  destination,  since  they  are  not  co-located.  In  this  section 
we  show  how  to  distribute  a  given  total  power  optimally  between  source  and  relay.  We  provide  first  a 
closed  form  analysis  in  the  ideal  case  where  there  are  no  decision  errors  at  the  relay  and  then  we  will 
show  some  performance  results  concerning  the  real  case  where  the  errors  are  taken  into  account. 
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1)  Error-free  S2R  link:  Under  the  assumption  that  there  are  no  errors  at  the  relay  side,  using  the 
same  derivations  introduced  in  Sec.  Ill,  the  optimal  detector  is  a  symhol-hy-symhol  detector  and  the 

signal-to-noise  ratio  on  the  /c-th  symbol  in  the  n-th  block  is 

42 

SNRk{n)  =  (a\Esd{k,k)\^  +  (1  -  Oi)\Kd{k,k)\^)  ,  (48) 

fox  k  =  ,N  and  n  =  i,  i  +  1.  The  error  probability  for  binary  antipodal  constellation,  conditioned 
to  a  given  channel  realization,  is  given  by 

Pe\h{k)  =  ^erfc  (y/0.5  5A^i?fc)  (49) 


where  SNRk  is  given  by  (48).  For  each  sub-carrier  k,  the  signal-to-noise  ratio  SNRk  is  given  by  the 
sum  of  two  statistically  independent  random  variables,  each  one  distributed  according  to  a  pdf  with 
two  degrees  of  freedom.  Thus,  using  (48)  and  (49),  the  BER  P),  averaged  over  the  channel  realizations 
is  given  by  [33] 


Pb 


1  72 

2  72  -  71 


IP) 


71 

7i  -  72 


(50) 


where 


A2  ^.^2 


aa 


71  :  = 


2  ^cx.  ’  *  2 


)  —  ^hiPh  +  !)• 


From  (50),  we  infer  that,  if  the  errors  in  the  relay’s  detection  are  negligible,  D-OSTC  scheme  achieves, 
as  expected,  the  maximum  available  diversity  gain,  equal  to  two.  The  optimum  value  of  a  can  be  found 
by  minimizing  (50).  Since  the  average  BER  (50)  is  a  convex  function  with  respect  to  a  [40],  the  optimal 
value  of  the  minimization  admits  a  unique  solution. 

It  is  straightforward  to  show  that,  if  D  is  equipped  with  ur  antennas,  the  achieved  diversity  gain  is 


2  Hr. 


2)  S2R  link  with  errors:  When  the  errors  at  the  relay  side  are  explicitly  taken  into  account,  it  is  not 
easy  to  derive  the  performance  of  the  optimal  detector  (46)  in  closed  form.  In  this  case,  it  is  interesting 
to  check  the  performance  of  the  sub-optimal  detector  (47),  to  quantify  the  loss  with  respect  to  the  more 
complex  but  optimal  detector  (46).  In  [40]  a  closed  form  expression  for  the  BER  of  the  sub-optimum 
scheme,  in  the  presence  of  relay  decision  errors,  was  derived.  For  a  given  total  transmitted  power  from 
S  and  R,  we  can  optimize  the  power  allocation  between  S  and  R,  depending  on  the  relative  distances 
between  S,  R  and  D,  in  order  to  minimize  the  final  average  BER.  We  address  this  issue  in  the  following 
example. 
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Example  1  -  Optimal  power  allocation:  We  show  now  the  behavior  of  the  final  average  BER  as  a  function 
of  the  power  allocation  between  S  and  R,  depending  on  the  relative  distances  between  S,  R  and  D.  As 
an  example,  in  Fig.  9  we  report  the  average  BER  vs.  a,  as  defined  in  Section  III-C.2,  for  different  values 
of  the  distance  drd  (and  thus  of  SNRn)  between  R  and  D  (all  distances  are  normalized  with  respect  to 
the  distance  d^d  between  S  and  D).  In  Fig.  9a),  we  consider  the  ideal  case  where  there  are  no  errors  at 
the  relay  node.  The  SNRd  at  the  final  destination  is  fixed  equal  to  10  dB.  We  can  observe  that,  when 
dsd  =  drd.  =  the  value  of  a  that  minimizes  the  average  BER  is  a  =  0.5,  i.e.  the  two  transmitter  use  the 
same  power.  However,  as  R  gets  closer  to  D,  the  optimal  a  tends  to  increase,  i.e.  the  system  allocates 
more  power  to  S,  with  respect  to  R.  The  reverse  happens  when  drd  is  greater  then  1.  Thus,  as  expected, 
the  system  tends  to,  somehow,  put  S  and  R  in  the  same  conditions  with  respect  to  D,  in  order  to  get 
the  maximum  diversity  gain. 

The  real  case,  where  there  are  decision  errors  at  the  relay  node,  is  reported,  as  an  example,  in  Fig.  9b), 
where  the  average  BER  is  again  plotted  as  a  function  of  a,  but  for  different  values  of  the  SNRr  at 
the  relay  node.  Interestingly,  we  can  observe  that,  as  SNRr  decreases,  the  system  tends  to  allocate  less 
power  to  the  relay  node  (the  optimal  value  of  a  is  greater  than  0.5),  as  the  relay  node  becomes  less  and 
less  reliable. 

D.  Performance 

In  this  section,  we  compare  alternative  cooperative  strategies.  We  assume  a  block  length  N  =  2,2  and 
channel  order  L  =  6.  To  make  a  fair  comparison  of  the  alternative  transmission  schemes,  we  enforce  all 
systems  to  transmit  with  the  same  overall  power.  More  specifically,  if  V  is  the  total  power  radiated  by 
the  non  cooperative  scheme,  we  denote  by  Vi  the  power  radiated  by  S  during  the  first  time-slot  and  by 
aVu  and  (1  —  a)Vii  the  power  spent  respectively  by  S  and  R  in  the  second  time-slot.  Since  the  overall 
radiated  power  is  always  V,  it  must  he.  V  =  Vj  +  Vn.  The  coefficient  a  is  chosen  in  order  to  minimize 
the  final  average  bit-error  probability  (50)  (see  also  Example  1)^^  The  power  Vi  is  chosen  in  order  to 
achieve  a  required  average  SNRr  at  the  relay,  defined  as  SNRr  :=  ■  All  distances  in  the  network 

are  normalized  with  respect  to  the  distance  dsd  between  S  and  D. 


'*We  use  theoretical  derivations,  valid  in  the  absence  of  errors  at  the  relay,  to  simplify  the  strategy.  One  could  improve  upon 
this  choice  by  using  the  BER  resulting  in  the  presence  of  errors  at  the  relay. 
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Fig.  9.  Average  BER  vs.  a;  a)  no  errors  at  R;  b)  including  errors  at  R. 


Example  2  -  ML  V5.  sub-optimum  detector:  In  Figs.  10  and  11,  we  compare  the  average  BER  obtained 
using  alternative  cooperative  and  non-cooperative  schemes.  The  BER  is  averaged  over  2000  independent 
channel  realizations.  All  curves  are  plotted  versus  the  SNR  in  D,  defined  as  SNRd  :=  ■  This  is 

also  the  SNR  of  the  single-hop  (non-cooperative)  case.  The  variance  of  the  noise  at  both  R  and  D  is 
unitary.  In  this  example,  we  set  dgr  =  0.1  and  d^d  =  0.9.  The  results  shown  in  Fig.  10  and  11  are  achieved 
transmitting  with  a  power  Vi  yielding  an  average  SNRr  at  the  relay  equal  to  15  dB,  for  all  values  of 
SNRd  reported  in  the  abscissas.  Since  the  noise  power  and  SNRr  are  both  fixed,  increasing  SNRd 
means  that  "P//  increases.  In  Fig.  10  and  11  we  report,  for  the  sake  of  comparison,  the  average  BER 
obtained  with  the  following  schemes:  i)  the  single -hop  method  (dotted  line);  ii)  the  ideal  ML  detector 
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for  0-DSTC  scheme,  with  no  errors  at  the  relay  (dashed  and  dotted  line);  iii)  the  real  ML  detector, 
incorporating  the  decision  errors  at  the  relay  (dashed  line);  iv)  the  suh-optimum  scalar  decoder  for  O- 
DSTC  scheme,  showing  both  the  theoretical  average  BER  (solid  line)  and  the  corresponding  simulation 
results  (circles),  obtained  using  a  Zero-Forcing  detector  at  the  relay  node.  In  Fig.  11  we  report  the 
performance  of  the  optimal  versus  sub-optimal  detector,  on  real  wideband  channel  measurements,  kindly 
provided  by  the  group  of  the  University  of  Bristol  led  by  Proff.  A.  Nix  and  M.  Beach.  The  data  are 
collected  in  a  urban  environment  (the  city  of  Bristol)  in  the  band  between  1920  and  1930  MHz. 


Fig.  10.  Comparison  between  average  BER  vs.  SNRd  (dB)  achieved  with  different  decision  schemes  over  Rayleigh  channels: 
Single  S-D  link  (dotted  line);  ideal  ML  (dashed-dotted  line);  real  ML  detector  (dashed  line);  sub-optimum  receiver  -  theoretical 
results  (solid  line)  and  simulation  (circles);  SNRr  =  15  dB. 

We  can  observe  a  very  good  agreement  between  our  theoretical  derivations  for  the  sub-optimum  detector 
and  the  corresponding  simulation  results.  The  floor  on  the  BER  of  the  sub-optimum  receiver  is  due  to  the 
decision  errors  at  the  relay  node.  It  is  also  interesting  to  notice,  from  Fig.  10,  that  the  sub-optimum  O- 
DSTC  scheme  exhibits  performance  very  close  to  the  optimal  0-DSTC  ME  detector,  at  low  SNRd,  i.e. 
before  the  BER  floor,  when  the  relay  is  relatively  close  to  the  source.  This  indicates  that  the  sub-optimum 
detector  is  indeed  a  very  good  choice,  under  such  a  scenario,  because  it  is  certainly  less  complicated  to 
implement  than  the  ME  detector.  Most  important,  differently  from  ME,  the  sub-optimum  scheme  does 
not  require  any  exchange  of  information  between  R  and  D,  about  the  BER  in  R.  The  price  paid  for  this 
simplicity  is  that  the  R  node  must  have  a  sufficiently  high  SNR  to  guarantee  that  the  BER  of  interest 
be  above  the  floor.  In  Figs.  10  we  have  also  reported  the  average  BER  (solid  line  with  starts)  obtained 
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Fig.  11.  Comparison  between  average  BER  vs.  SNRd  (dB)  achieved  with  different  decision  schemes  over  real  channel 
measurements:  Single  S-D  link  (dotted  line);  ideal  ML  (dashed-dotted  line);  real  ML  detector  (dashed  line);  sub-optimum 
receiver  -  theoretical  results  (solid  line)  and  simulation  (circles);  SNRh  =  15  dB. 


using  a  transmission  strategy,  for  0-DSTC  scheme,  where  instead  of  OFDM,  in  the  S/R  slot  we  used 
a  linear  precoding  method  that  insures  minimum  BER  at  the  relay,  under  the  assumption  of  adopting  a 
(suhoptimal)  MMSE  linear  decoder  (solid  line  with  stars).  In  such  a  case,  we  observe  that,  with  minimum 
additional  complexity  at  the  relay,  the  performance  of  the  suhoptimal  0-DSTC  scheme  becomes  closer 
to  the  ME  decoder  because  of  the  lower  BER  at  the  relay. 

Einally,  looking  at  the  slopes  of  the  average  BER  curves  of  the  ME  0-DSTC  detector,  shown  in  Eigs.  10 
and  1 1,  it  is  worth  noticing  that,  in  the  absence  of  errors  at  the  relay,  the  cooperative  scheme  achieves  full 
spatial  diversity  gain,  provided  that  the  relay  can  be  used,  as  we  have  assumed  in  this  section.  In  practice, 
there  are  two  reasons  for  the  lack  of  full  diversity  gain.  The  first  one  is  that,  in  a  network  where  the 
relays  are  randomly  spatially  distributed,  the  probability  of  finding  no  relay  is  nol  zero.  The  second  one  is 
fhat,  given  a  sef  of  ferminals  available  as  relays,  because  of  fhe  presence  of  decoding  errors  af  fhe  relay 
side,  all  cooperative  schemes  exhibit  an  asymptotic  average  BER  behavior  proportional  to  1/SNRd- 
Nevertheless,  there  is  a  considerable  coding  gain,  which  justifies  fhe  use  of  cooperation.  Indeed,  a  more 
attentive  look  at  the  results  shows  that  the  average  BER  starts  approaching  the  slope  with  maximum 
diversity,  as  far  as  the  errors  at  the  relay  are  negligible  with  respect  to  the  errors  at  the  destination.  Then, 
when  the  errors  at  the  relay  become  dominant,  the  final  BER  curve  follows  fhe  1/ SNR  behavior. 
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E.  Choice  of  the  Constellation  Order  in  the  Source-Relay  slot 

The  other  major  critieal  aspect  of  cooperative  schemes  is  their  rate  loss  due  to  the  insertion  of  the 
S/R  time  slot.  As  an  example,  if  all  the  links  would  use  a  BPSK  constellation,  the  rate  loss  factor 
would  he  1/2.  To  reduce  this  loss  factor,  we  can  use  higher  order  constellations  in  the  S/R  link,  with 
respect  to  the  constellations  used  in  the  other  links,  so  that  the  duration  of  the  S/R  slot  can  he  made 
smaller  than  the  duration  of  the  other  slots.  In  this  section,  we  assume  BPSK  transmissions  over  all 
links,  except  the  S/R  link,  where  the  constellation  order  is  allowed  to  increase.  More  specifically,  using 
a  constellation  A  of  cardinality  M  =  2^'’  in  the  S/R  link,  the  rate  loss  factor  is  rib /{rib  +  1)-  On  the 
other  hand,  cooperation  reduces  the  final  BER  and  then  it  induces  a  capacity  increase.  To  quantify  the 
overall  balance  in  terms  of  rate,  we  compared  the  maximum  rate  achievable  by  0-DSTC  system  with 
the  maximum  rate  achievable  with  a  non-cooperative  scheme.  We  define  as  achievable  rate  the  maximum 
number  of  bits  per  symbol  (bps)  that  can  be  decoded  with  an  arbitrarily  low  error  probability,  provided 
that  sufficient  error  correction  coding  is  incorporated  in  the  system,  conditioned  to  the  assumptions  (al)- 
(a6)^^.  We  have  shown  in  Section  III  that  the  combination  of  0-DSTC  and  OFDM  makes  the  overall 
time-dispersive  channel  equivalent  to  a  set  of  parallel  non-dispersive  sub-channels.  The  final  S2D  link, 
over  each  sub-channel,  in  the  presence  as  well  as  in  the  absence  of  the  relay  link,  can  always  be  made 
equivalent  to  a  binary  symmetric  channel  (BSC)  with  cross-over  probability  depending  on  the  specific 
cooperative  (or  non-cooperative)  scheme  adopted.  Thus,  the  maximum  rate  R{k\h)  that  can  be  reliably 
transmitted,  over  the  /c-th  sub-carrier,  for  a  given  channel  realization  h,  incorporating  the  rate  loss  due 
to  the  insertion  of  the  S2R  slot,  is 

^(^1^)  =  ^  x/nfc  *^BSc(Pe|h(k))  bps.  (51) 

where  Pe\h{k)  denotes  the  binary  error  probability  on  the  /c-th  sub-carrier,  conditioned  to  the  channel 
realizations,  Cbsc{p)  '■=  1  +plog2(p)  -|-  (1  —  p)log2(l  —  p)  :=  1  —  H{p)  is  the  capacity  of  a  binary 
symmetric  channeP^  with  crossover  probability  p.  From  (51)  we  infer  that,  because  of  the  S2R  link, 
cooperative  transmission  induces  a  systematic  rate  loss  of  rib /{rib  +  1)>  with  respect  to  the  case  of  no 
cooperation.  But,  at  the  same  time,  cooperation  yields  a  smaller  error  probability  Pe\h{k)  and  thus  a 
higher  CBSc{Pe\h{k))-  Then,  we  may  expect  a  trade-off  in  the  choice  of  rib-  This  trade-off  can  be  better 
understood  through  the  following  example. 

'^It  is  important  to  remark  that  the  rate  defined  above  is  smaller  than  the  capacity  of  the  system,  because  the  proposed  scheme 
is  designed  to  maximize  the  spatial  diversity  gain  and  not  to  maximize  information  rate. 

*^We  can  use  this  formula  because  the  S-D  is  always  BPSK,  regardless  of  the  constellation  used  in  the  S2R  link. 
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Fig.  12.  Achievable  rate  (bps)  vs.  SNR  (dB)  -  non  cooperative  case  (solid  line),  cooperative  case  using:  BPSK  (circle  marker), 
QPSK  (star  marker),  16-QAM  (square  marker),  64-QAM  (‘+’),  in  the  S2R  link  -  a)  SNRr  =  15  dB;  b)  SNRr  =  3  dB. 


Example  3.  Rate  and  diversity  gain:  As  an  example,  we  report  in  Fig.  12  a)  and  b)  the  achievable  rate 
vs.  the  SNRd  in  D,  for  an  SNRr  in  R  equal  to  15  and  3  dB,  respectively,  for  different  choices  of 
the  constellation  used  in  the  S2R  link,  achieved  with  or  without  cooperation.  To  preserve  the  receiver 
simplicity,  zero-forcing  equalization  and  symbol-by-symbol  detection  are  performed  at  the  relay.  We 
can  see  that,  at  high  SNRr,  the  non-cooperative  case  approaches  the  maximum  value,  equal  to  1  bps, 
whereas  the  cooperative  cases  tend  to  an  asymptote  less  than  1,  depending  on  the  constellation  used  in 
S2R  slot.  We  observe  from  Fig.l2  a)  that,  for  SNRr  =  15  dB,  increasing  the  constellation  order  from 
BPSK  to  16-QAM  in  the  S2R  link  improves  the  achievable  rate;  however  passing  from  16  —  QAM 
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to  64  —  QAM  does  not  induce  any  further  gain  because  of  the  higher  BER  at  the  relay.  For  lower 
values  of  SNRr,  i.e.  SNRr  =  3  dB  for  example,  there  is  no  appreciable  rate  gain  in  increasing  the 
constellation  order  because  of  the  excessive  BER  at  the  relay.  Nevertheless,  it  is  interesting  to  notice  that, 
at  low/medium  SNRd  at  the  final  destination  (within  a  range  depending  on  SNRr),  the  cooperative 
case  can  outperform  the  non-cooperative  case  also  in  terms  of  achievable  rate,  because  the  BER  decrease 
can  more  than  compensate  the  rate  insertion  loss  due  to  the  S2R  slot. 

E  Synchronization 

Besides  errors  at  the  relay  node,  one  more  distinguishing  feature  of  DSTC  is  that  the  cooperating 
transmit  antennas  are  not  co-located.  This  means  that  the  packets  arriving  at  the  final  destination  from 
source  and  relays  might  be  asynchronous.  Interestingly,  if  the  difference  in  arrival  times  is  incorporated 
in  the  CP  used  from  both  S  and  R,  D  is  still  able  to  get  N  samples  from  each  received  block,  without 
inter-block  interference  (IBI).  In  such  a  case,  the  different  arrival  time  does  not  cause  any  trouble  to 
the  final  receiver.  In  fact,  let  us  take  as  a  reference  time  the  instant  when  the  i-th  block  coming  from 
R  arrives  at  D.  If  the  block  coming  from  S  arrives  with  a  delay  of  samples,  the  only  difference 
with  respect  to  the  case  of  perfect  synchronization  is  that  the  transfer  function  Asd{k)  in  (44)  will  be 
substituted  by  .  From  (44),  it  is  clear  that  such  a  substitution  does  not  affect  the  useful 

term,  as  it  only  affects  the  interfering  term.  However,  in  the  hypothesis  of  Rayleigh  fading  channel, 
Ksd{k)  is  statistically  indistinguishable  from  Ksd{k)e~^‘^'"^'^^/^ .  Hence,  the  combination  of  Alamouti 
(more  generally,  orthogonal  STC)  and  OFDM  is  robust  with  respect  to  lack  of  synchronization  between 
the  time  of  arrival  of  packets  from  S  and  from  R  (as  long  as  (a6)  holds  true).  The  price  paid  for  this 
robustness  is  the  increase  of  the  CP  length  L,  which,  in  its  turn,  reflects  into  a  rate  loss.  However,  this 
loss  can  be  made  small  by  choosing  a  blocklength  N  much  greater  than  L  or  by  selecting  only  relays 
that  are  relatively  close  to  the  source,  so  as  to  make  the  relative  delay  small. 

IV.  Regenerative  vs.  non-regenerative  relays  in  single-user  scenario 

In  non-regenerative  schemes,  the  relay  node  simply  amplifies  and  retransmits  the  received  signal,  with¬ 
out  performing  any  A/D  conversion  on  the  signal.  Thus,  A&F  can  be  useful  to  simplify  the  implementation 
of  the  relay,  because  the  relay  of  an  A&F  system  only  needs  to  have  an  antenna  and  a  RF  amplifier. 

Since  no  detection  can  be  performed  at  the  relay  side  in  the  A&F  scheme,  the  relay  can  only  retransmit 
the  received  signal.  Thus,  in  order  to  implement  the  distributed  version  of  Alamouti  scheme,  in  the  first 
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time-slot,  S  has  to  transmit,  consecutively,  the  blocks  — and  AsGs*{iY'^.  In  the  second  time- 
slot,  S  transmits,  consecutively,  AsFs{i)  and  then  ylsGs*(z-|-l)  and  R  sends  Ar{— H srF s{i+l)  +v ^{i)) 
first  and  then  Ar{HsrGs*{i)  +  vji{i  +  1)).  The  amplitude  coefficients  As  and  Ar  are  used  to  impose 
the  power  available  at  the  S  and  R  nodes,  respectively.  In  the  A&F  case,  differently  from  the  D&F  case, 
the  coefficient  Aj.  depends  also  on  the  S2R  channel  as  well  as  on  the  noise  at  the  R  node.  Clearly,  A^ 
changes  depending  on  which  strategy  is  implemented  in  R. 

Thanks  to  the  combination  of  Alamouti’s  coding  and  OFDM,  the  overall  systems  is  equivalent  to  a 
series  of  N  parallel  channels.  Proceeding  as  in  Sec.  Ill  to  obtain  (44),  the  received  symbol  pertaining 
to  the  k-th  sub-carrier,  in  the  i-th  and  i  +  1-th  time  slots  is  given  by^^ 


rk  = 


9k  0 

0  gk 


^k  T  ^ki 


(52) 


where  Sk  =  [sk{i),Sk{i  +  1)]'^,  Qk  :=  Al\Asd{k)Y  +  Al\Ard{k)Y\-l^sr{k)Y  and  Uk  is  a  Gaussian 
vector  with  zero  mean  and  diagonal  covariance  matrix  C„  =  all,  with  a‘1  =  {A^\Ky.d{k)Y a^  +  a^) 

{Al\Ksd{k)?  +  Al\Kd{k)\Msr{k)Y)- 


Example  3.  Comparison  between  Ak.F  and  D&iF:  In  Fig.  13,  we  compare  the  average  BER  vs.  the 
SNRd  at  the  destination  node,  obtained  using  the  following  strategies:  a)  Decode  and  forward  using 
ML  detector  (dashed  line)  or  sub-optimal  detector:  Theoretical  value  (solid  line)  and  simulation  results 
(circles);  b)  amplify  and  forward  (dashed-dotted  line);  c)  single  hop  (non-cooperative)  case  (dotted  line). 
The  block  length  is  =  16;  the  channels  are  simulated  as  FIR  filters  of  order  Lh  =  6,  whose  taps  are 
iid  complex  Gaussian  random  variables  with  zero  mean  and  variance  1/d^.  The  SNRji  at  the  relay  is 
equal  to  20  dB.  Comparing  the  D&F  and  A&F  schemes,  we  observe  that  the  D&F  method  performs 
better  than  the  A&F  at  low  and  intermediate  SNRd  values,  but  for  high  values  of  SNRd,  the  A&F 
performs  better.  This  shows  that  A&F  is  indeed  a  valuable  choice. 

V.  Comparison  among  alternative  STC  techniques  in  a  multi-user  context 

In  this  section,  we  compare  alternative  DSTC  strategies  in  a  multi-user  scenario.  We  consider  a  cell 
of  radius  300  m  with  Ntot  =  200  total  radio  terminals,  located  randomly.  We  consider  only  the  uplink 

'"^Differently  from  A&F  scheme,  in  the  D&F  system  there  is  no  constraint  on  the  sequence  of  information  blocks  transmitted 
from  the  source  node. 

'^We  drop  the  block  index  i  for  simplicity  of  notation,  because  the  same  relationships  hold  true  for  all  blocks. 
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Fig.  13.  Average  BER  vs.  SNRd  (dB)  achieved  with  different  strategies:  a)  Decode  and  forward  using  ML  detector  (dashed 
line)  or  sub-optimal  detector:  Theoretical  value  (solid  line)  and  simulation  results  (circles);  b)  amplify  and  forward  (dashed-dot 
line);  c)  single  hop  (non-cooperative)  case  (dot  line);  SNRr  =  20  dB. 


channel  and  we  assume  that  the  base  station  (destination)  has  two  real  antennas.  Within  the  set  of  all 
radio  terminals,  =  10  is  the  number  of  sources,  whereas  the  remaining  nodes  are  potential  relays. 
All  channels  are  slowly-varying,  Rayleigh  flat  fading.  The  alternative  strategies  are  compared  enforcing 
the  same  overall  radiated  energy  and  the  same  bit  rate  in  the  SR2D  phase.  In  case  of  cooperation,  the 
energy  includes  the  energy  used  to  send  data  from  the  source  to  the  relays  and  the  sum  of  the  energy 
used  by  source  and  relay  to  transmit  their  data  to  the  destination.  For  each  channel  realization  and  radio 
nodes  distribution,  we  associate  a  relay  to  a  source  according  to  the  protocol  described  in  Section  II-D. 
We  have  used  a  16-QAM  constellation  for  the  conventional  SISO  system,  that  acts  as  a  benchmark  term, 
while  for  the  DSTC  schemes  the  following  choices  have  been  made,  in  order  to  enforce  the  same  bit 
rate  in  the  SR2D  phase: 

a)  D-OSTC.  As  the  symbol  rate  of  the  OSTC  scheme  is  the  same  as  a  SISO  system,  the  constel¬ 
lation  used  is  the  same,  both  if  the  sources  finds  a  relay,  and  if  it  transmits  alone.  We  have 
chosen  16-QAM  so  that  each  source  transmits  8  bits  every  22  symbol  intervals.  Two  symbol 
intervals  are  reserved  for  the  common  S2R  phase  {Ts2R  =  2)  hence  also  in  this  phase  16-QAM 
is  adopted. 

b)  D-BLAST.  A  4-QAM  constellation  is  used  in  case  of  cooperation  and  a  16-QAM  constellation 
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is  used  in  case  of  no-cooperation. 

The  cooperative  case  can  use  a  lower  order  constellation  because  a  2  x  2  D-BLAST  provides  a  higher 
transmission  rate. 

c)  D-FDFR.  A  4-QAM  constellation  is  used  in  case  of  cooperation  and  a  16-QAM  constellation 
is  used  in  case  of  no-cooperation. 

Also  in  this  case,  as  with  D-BLAST,  the  cooperative  case  can  use  a  lower  order  constellation  with  respect 
to  the  non-cooperative  case. 


V-BLAST  4-QAM,  SNIR  at  the  Relay  =12.5  dB 
-O-  Full  Rate-Full  Div  4-QAM,  SNIR  at  the  Relay  =1 2.5  dB  | 

-  No  Coop  1 6-QAM 

-0-  Orthogonal  STC  1 6-QAM,  SINR  at  the  Relay  =1 2.5  dB  I 
-A-  Full  Rate-Full  Div  4-QAM,  SINR  at  the  Relay  =1 5  dB 
Orthogonal  STC  1 6-QAM,  SINR  at  the  Relay  =1 5  dB 
V-BLAST  4-QAM,  SNIR  at  the  Relay  =15  dB 


Averages  over: 

5  realizations  of  the  Terminals’  deployment 
1 0(  X  5  )  realizations  of  the  Active  Sources  locations 
20(  X  10  X  5 )  channels  realizations  (for  each  Source) 
Block  Length :  Nb  =  1 0 

Number  of  Terminals  =200 
Number  of  Active  Sources  =  1 0 
a  =  0.1 


5  10  15  20  25 

Average  SNR  at  the  Destination 


Fig.  14.  BER  comparison  of  different  DSTC  schemes 

The  simulation  results  are  reported  in  Fig.  14.  Throughout  the  simulations,  the  power  used  in  the  S2R 
slot  is  one  tenth  of  the  power  used  in  the  non-cooperative  case.  As  explained  in  the  last  section,  this 
portion  could  be  optimized.  A  node  is  chosen  as  a  relay  if  its  Signal  to  Noise  plus  Interference  Ratio 
(SNIR),  conditioned  to  the  channel,  exceeds  a  threshold  of  12.5  or  15  dB.  Requiring  a  SNIR  of  12.5 
dB,  a  relay  has  been  found  with  probability  pr{1)  =  0.72,  whereas  for  SNIR  =  15  dB,  we  obtained 
Pr{1)  =  0.65.  Clearly,  increasing  the  target  SNIR,  it  decreases  the  probability  of  finding  a  relay,  but,  at 
the  same  time  there  are  less  decision  errors  at  the  relay.  The  overall  performance  is  then  a  combination 
of  these  two  aspects.  The  average  BER  reported  in  Fig.  14  takes  into  account  both  situations  where  the 
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relay  has  been  found  or  not.  We  can  check  from  Fig.  14  that  indeed,  increasing  the  SNIR  from  12.5  to 
15,  even  though  Pro(l)  decreases,  the  floor  on  the  BER  decreases  by  more  than  a  decade.  Of  course, 
this  result  is  also  a  consequence  of  the  relay  density.  Finally,  in  the  case  of  a  SNIR  of  15  dB,  we  can 
observe  a  gain  of  approximately  3  dB  at  BER  =  10“^. 

VI.  Cooperative  scheme  composed  of  the  cascade  of  SISO  and  virtual  MIMO  links 

In  this  section,  we  focus  on  systems  where  the  final  destination  has  multiple  antennas.  This  means  that, 
if  more  relays  cooperate  with  the  source  in  sending  information  to  the  destination,  we  may  establish  a 
virtual  MIMO  link  between  the  relays  and  the  destination.  The  specific  goal  of  this  section  is  to  derive  the 
energy  allocation  and  time  sharing  that  provide  a  performance  gain  with  respect  to  the  non-cooperative 
case. 

A  non  negligible  issue  is  the  geometry  of  the  system,  i.e.,  the  relative  location  of  the  nodes  involved  in 
the  communication.  Aiming  at  establishing  achievable  rate  limits,  this  aspect  has  been  recently  explicitly 
taken  into  account,  in  [70]  for  the  Gaussian  channel,  in  [71],  for  fading  channels  in  the  low-power  regime, 
and  in  [72],  where  both  an  ergodic  capacity  analysis  and  an  outage  probability  analysis  are  carried  out 
for  a  one-relay  channel,  considering  the  effect,  on  the  performance  limits,  of  the  relay  location. 

We  consider  a  Rayleigh  flat-fading  channel  in  a  scenario  in  which  a  source  terminal  exploits  >  1 
relays  to  form  a  virtual  array,  thus  enabling  a  MIMO  transmission  towards  a  final  destination,  which  is 
equipped  with  >  1  antennas.  A  sketch  of  the  considered  geometry  is  depicted  in  Fig.  15. 


Fig.  15.  Transmission  through  a  virtual  array  of  relays. 

In  this  section,  we  establish  the  theoretical  performance  limits,  both  in  the  ergodic  and  non-ergodic 
case,  emphasizing  the  dependence  of  these  limits  on  the  following  parameters:  i)  the  distances  between 
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source,  relays,  and  destination;  ii)  the  power  allocation  and  iii)  the  time  allocation,  between  the  two 
phases. 

In  the  general  case,  the  destination  can  receive  signals  in  both  the  transmission  phases.  Furthermore, 
the  source  may  partecipate  to  the  MIMO  phase,  thus  contributing  to  the  system  diversity.  For  simplicity, 
however,  we  will  assume  that  the  source  transmits  in  the  first  phase  only,  and  the  destination  receives 
in  the  second  phase  only.  Thus,  we  have  the  cascade  of  a  one-to-many  link  between  the  source  and  the 
relays,  which  we  will  call  a  n^-SISO  channel,  and  a  MIMO  link  between  the  virtual  array  of  relays 
and  the  physical  array  at  destination.  The  significance  of  our  resulf  is  nof  affecfed  by  fhis  simplifying 
assumption,  since  if  clearly  represents  a  worst-case  scenario,  with  respect  to  the  general  case. 

Throughout  this  section,  we  shall  assume  decode  &  forward,  [53],  as  the  relaying  mode.  Our  goal  is 
to  show  that  a  MIMO  communication  can  be  effectively  implemented  with  single-antenna  transmitters, 
thus  yielding  the  well  known  benefits  of  MIMO  systems.  Special  emphasis  will  be  given  to  the  role  of 
the  geographical  distribution  of  the  relays,  power  allocation  and  time-sharing. 

Both  in  the  ergodic  and  non-ergodic  cases,  we  will  compare  the  performance  limits  of  the  proposed 
system  with  those  of  a  traditional  direct  source-destination  SIMO  link. 

1 )  System  model  and  problem  formulation:  We  assume  that  the  source  is  equipped  with  a  single 
physical  antenna,  whereas  the  destination  is  equipped  with  an  array  of  n^)  antennas.  In  between,  there  is 
a  set  of  Ur  single-antenna  relays,  that  assist  the  source  in  the  transmission,  through  the  formation  of  a 
virtual  array.  In  our  setup  the  relays  do  not  have  their  own  data  to  send.  The  communication  is  divided  in 
two  phases:  in  the  first  phase,  that  we  call  the  “source  to  relays”  (S2R)  phase,  the  source  sends  its  data 
to  the  relays  in  a  broadcast  fashion.  The  relays  receive  and  decode  the  data  and,  in  the  second  phase, 
called  the  “relays  to  destination”  (R2D)  phase,  transmit  the  data  to  the  destination  as  in  a  MIMO  link. 
Since  the  scope  of  our  analysis,  at  this  stage,  is  not  to  find  the  best  suitable  transmission  strategy,  but 
to  investigate  the  information-theoretical  limits  of  the  proposed  protocol,  we  do  not  assume  a  specific 
space-time  code  for  the  R2D  phase.  However,  we  will  assume  that  in  order  to  achieve  full  diversity,  all 
the  relays  must  decode  the  entire  set  of  data  in  the  first  phase,  thus  enabling  the  use  of  full  diversity 
schemes  in  the  second  phase. 


We  assume  a  Rayleigh  flat-fading  channel  in  all  the  links.  The  discrete-time  channel  between  the  source 
and  one  of  the  relays  (S2R),  is  modeled  by  the  input/output  relation 


yr,k  — 


-hkX  +  W, 


where  Vs  is  the  transmission  power,  hk  is  a  complex  gaussian  random  variable  with  zero  mean  and  unit 
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variance  that  represents  the  normalized  channel  coefficient,  dk  is  the  distance^^  between  the  source  and 
the  k-th  relay,  tj  is  the  path  loss  exponent,  x  is  the  transmitted  symbol,  assumed  to  have  unit  energy, 
i.e.  =  1,  and  finally  w  is  an  addifive  while  complex  Gaussian  noise  term  wilh  zero  mean  and 


variance  cj^. 

In  fhe  R2D  phase,  Ihe  MIMO  link  between  the  relays  and  the  destination  is  modeled  by  the  input/output 
relation 


z  = 


HDy  +  w, 


(54) 


where  z  is  the  n^-size  received  vector,  y  is  a  vector  whose  k-th  entry  represents  the  symbol  transmitted 
by  the  k-th  relay,  which  depends  on  what  the  relay  has  received  in  the  first  phase,  Vr  is  the  transmit 
power,  which  is  equally  divided  among  the  relays,  w  is  an  UD-size  additive  noise  vector  with  independent 
zero  mean  circularly  symmetric  Gaussian  entries  of  variance  cj^,  H  is  the  ur  x  ur  normalized  channel 
matrix  ,  whose  elements  are  i.i.d.  circularly  symmetric  complex  Gaussian  random  variables  with  zero 
mean  and  unit  variance,  finally,  D  is  a  ur  x  ur  diagonal  malrix  lhat  takes  into  account  the  average 
path  losses,  and  is  defined  as 

Ddiag  ^(^1 -\- )  •  •  • ;  (l  + (^^) 


where  k  =  1, . . .  ,nR,  are  fhe  dislance  between  the  relays  and  the  destination,  and  t]  is  the  path 

loss  exponent,  usually  in  the  interval  [2,6]. 

We  assume  that  the  source  sends  a  total  of  Mi,  bits  to  the  destination  in  an  interval  of  T  channel  uses, 
i.e.,  at  an  average  transmission  rate  of  i?  =  M^/T  bits  per  channel  use  (bpcu).  Ts2R  channel  uses  will 
be  employed  for  the  S2R  phase,  whereas  Tr2d  channel  uses  are  reserved  for  transmission  between  the 
relays  and  the  destination.  We  define  e  as  fhe  fofal  energy,  comprehensive  of  bofh  fhe  S2R  and  fhe  R2D 
phases,  employed  fo  send  fhe  data  to  destination.  The  average  power  is  given  by 

Vav^.  (56) 


We  assume  that  a  fraction  a  of  the  total  energy  is  employed  in  the  first  phase,  and  the  remaining  fraction 
(1  —  a)  is  reserved  for  the  second  phase.  The  power  term  Vs  in  (53)  becomes 


Vs  =  a 


Ts2R 


=  a 


Ts2rJ 


Va 


(57) 


'^Notice  that  with  the  propagation  model  assumed  in  (53),  the  distance  is  normalized  so  that  a  unit  distance  corresponds  to 
the  distance  at  which  the  transmission  power  is  halved. 
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Since  in  the  S2R  phase  the  relays  have  to  receive  Mb  hits,  the  rate  of  SISO  links  of  the  S2R  phases  is 


given  hy 


RS2R7^ -  =  7^ - R- 

4  S2R  4  S2R 


Similarly  to  (57),  we  define  the  average  power  employed  in  the  R2D  phase,  that  appears  in  (54),  as 


=  (1  -  a) 


Vd  =  {1-  a) 


whereas  the  rate  of  the  virtual  MIMO  channel  is 

RS2R7^  =  t^R.  (60) 

4  R2D  4  P12D 

In  the  following,  we  will  assess  the  performance  limits  of  this  communication  system,  comparing  them 
with  those  achievable  with  a  standard  direct  link  between  the  source  and  the  destination.  We  will  consider 
two  different  channel  evolution  models,  i.e.  ergodic  and  non-ergodic  fading  channels.  Transmission  over 
ergodic  channels  means  that  the  channel  coherence  time  is  sufficiently  short  to  encode  the  information 
in  codewords  that  last  across  many  independent  channel  realization.  Clearly,  the  channel  coherence  time 
should  in  any  case  be  sufficiently  long  to  enable  CSI  acquisition  at  the  receiver.  For  ergodic  channels, 
it  is  possible  to  define  the  capacity  in  the  classical  formulation,  as  the  maximum  information  rate  that 
the  channel  can  support  enabling,  in  principle,  an  almost  error-free  communication.  Non-ergodic  fading 
channels  refer  instead  to  a  very  slow  fading  environment  in  which,  due  to  delay  constraints,  coding  over 
multiple  channel  realization  is  impossible.  Thus,  for  every  channel  realization  it  is  possible  to  define  a 
supported  rate,  or  “instantaneous  capacity”,  but  it  is  known  that  the  capacity  of  the  channel  ,  in  the 
Shannon  sense,  is  zero.  In  this  case  the  performance  limits  are  conveniently  expressed  in  terms  of  outage 
probability,  with  reference  to  a  given  target  information  rate  R,  or  outage  capacity,  with  reference  to 
a  given  outage  probability  Pout-  The  first  is  defined  as  the  probability  that  the  instantaneous  capacity, 
which  is  a  random  variable,  is  less  than  the  required  rate.  The  second  is  the  rate  that  guarantees  that  the 
probability  that  a  single  channel  realization  can  support  it  is  less  or  equal  to  the  desired  Pout- 


The  scope  of  our  analysis  is  to  emphasize  the  role  played  by  the  parameters  introduced  in  (53)-(60), 
namely:  the  (relative)  distances  between  the  relays  and  the  source/destination,  d^,  dD,fc)  k  =  1,  ,nR,  the 
power  allocation  across  the  two  phases,  represented  by  the  parameter  a,  and  the  time  division  of  the  two 
phases  represented  by  the  ratio  T/Ts2R- 
We  examine  now  the  ergodic  case. 
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2)  Ergodic  capacity  analysis:  If  the  system  requirements  on  the  deeoding  delay  allow  for  a  coded 
hlocklength  which  is  much  greater  than  the  channel  coherence  time,  each  block  has  the  possibility  to  be 
reliably  decoded,  if  the  transmission  rate  is  below  the  ergodic  channel  capacity.  For  the  S2R  phase,  the 
capacity  of  the  link  between  the  source  and  the  k-th  relay,  measured  in  bits  per  channel  use,  is  given  by 

Cs2R,k  =  I  ’  k  =  1,. . .  ,nii,  (61) 

where  \hkf  /  (l  +  d^l)  >  k  =  1, . . .  ,nR,  are  the  path  losses  corresponding  to  a  single  channel  re¬ 
alization  on  the  riR  links.  It  is  easy  to  chek  that  the  path  losses  are  exponentially  distributed  random 
variables  with  parameters  (l  -|-  ,  k  =  1, ,  hr,  respectively. 

Since  for  an  end-to-end  reliable  communication  it  is  necessary  that  all  the  relays  perfectly  decode  the 
transmitted  codeword,  we  can  define  an  overall  capacity  in  the  S2R  link,  which  is  the  minimum  among 
all  the  capacities  of  the  ur  parallel  channels 

Cs2R  =  min  {Cs2R,l,  •  •  •  ,  Cs2R,nR}  ■  (62) 

In  the  D2R  phase,  we  have  a  MIMO  link  whose  ergodic  capacity  is  given  by 

CR2D  =  EH!^logdet(^Im  +  ^W^Y  (63) 

where  mmin  {ur,  no),  Im  is  the  Identity  matrix  of  size  m,  and  W  is  defined  as 

f  HDiHDf  ,  if  hr  >  no 

W{  ^  >  ,  J  R-  U  ^ 

[  {HDfHD,  if  UR  <  no 

We  have  that  if  Cs2R  >  Cr2d,  the  maximum  information  that  could  be,  in  principle,  transferred  to  the 
relays,  cannot  flow  toward  the  destination  because  the  MIMO  channel  is  not  able  to  support  it.  Viceversa, 
if  Cr2d  >  Cs2R’  it  is  clear  that  the  bottleneck  for  the  information  flow  lies  in  the  first  link.  The  capacity 
of  the  composite  link  we  are  considering  is  hence  given  by  the  minimum  among  the  two: 

C S2D  =  {C S2R,  C R2d}  ■  (65) 


For  the  sake  of  simplicity,  let  us  assume  that  the  relays  are  located  at  the  same  distance  from  the  source, 
i.e.  dk  =  dsR,  k  =  1, . . .  ,nR,  and  also  from  the  destination,  i.e.  dD,k  =  dRo,  k  =  I, . . .  ,nR. 

We  also  assume  that  the  relays  lay  on  the  direction  connecting  the  source  to  the  destination,  as  depicted 
in  Fig.  15,  i.e. 

dRD  +  dsR  =  dsD-  (66) 


In  this  case  we  have: 


VF  =  (1  +  {dsD  -  dsR)^)-^  W, 


(67) 
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where 


W 


HH 


H 


(68) 


if  UR  >  no 
H^H,  if  UR  <  no 

From  (61)  and  (63)  we  can  rewrite  the  capacities  in  the  two  phases,  hy  making  explicit  the  dependence 
on  the  distances  (Ird,  dsR,  and  dsn,  on  the  power  allocation  and  time-sharing,  represented  hy  a  and 
the  ratio  T/Ts2R,  respectively,  thus  obtaining 


^  TP 

^S2R  —  m 


log  I  1  +  ( 

J^S2RJ 


\h[ 


0-2  (l  + 


Cr2d  =  Eh  {  log 


J-ra  H“ 


((1  -  a) 


Va 


n-Rcr"^  (1  -I-  {dsD  -  dsR)^) 


W 


(69) 


(70) 


where  we  have  introduced  the  scale  factors  Ts2r/T  and  T£)2r/T ,  for  a  fair  comparison  with  a  direct 
transmission  system.  We  can  see  that,  for  a  given  distance  dsn  between  the  source  and  the  destination, 
the  effect  of  an  increase  in  dsR  is  to  increase  Cr2d  and  to  decrease  Cr2d,  thus,  we  can  expect  that 
there  will  be  a  value  d^^  of  the  distance  for  which  the  two  capacities  are  equal  and  yield  the  maximum 
achievable  rate  of  the  system.  Furhermore,  the  capacity  is  also  increased,or  decreased,  depending  on  the 
power  allocation  and  the  transmission  time  allocation. 


Fig.  16.  Ergodic  capacity  of  the  cooperative  system,  compared  with  the  capacity  of  a  direct  S— i-D  transmission  system 

In  Fig.  16  we  have  plotted  the  ergodic  capacity  of  a  cooperative  system  (blue  surface)  formed  by 
Ur  =  6  relays  and  a  receiver  with  =  3  antennas.  The  black  surface  represent  the  ergodic  capacity 
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of  the  direct  system,  the  red  curve  represents  the  distance  that  maximizes  the  capacity,  for  each  value  of 
a,  whereas  the  region  of  the  (a,  plane  comprised  between  the  two  hlack  contours  is  the  region  of 
distances  for  which  the  cooperative  system  outperforms  the  direct  system.  These  results  were  obtained 
with  dsD  =  100,  p  =  3,  and  Vavl^^^  =  100  dB ,  and  a  time-sharing  factor  Ts2r/T  =  0.5. 

We  can  see  that  there  exist  range  of  values  of  a  and  dsR  providing  a  significant  benefit  over  the 
SIMO  channel  between  source  and  destination,  thus  justifying  the  use  of  cooperative  relaying.  It  is 
important,  however  to  choose  properly  the  value  of  a  for  a  given  distance  (or  viceversa)  to  achieve  the 
best  performance.  Similar  results  would  be  obtained  for  a  fixed  a,  and  different  values  of  the  time-sharing 
factor. 

3)  Outage  based  analysis:  Depending  on  the  channel  coherence  time,  it  could  be  impractical  to  send 
code -blocks  whose  length  allows  to  achieve  the  ergodic  capacity  C  of  the  fading  channel.  In  a  slow 
fading  scenario,  it  is  frequent  to  consider  the  single  channel  realizations  and  evaluate  the  perfomance  in 
terms  of  outage  capacity  or  outage  probability.  In  practice,  each  channel  realization  yields  a  realization 
of  a  random  variable,  called  “instantaneous  capacity”,  Ch  (see  for  example  [59]),  which  is  the  capacity, 
in  the  Shannon  sense,  that  one  would  have  if  that  channel  realization  had  an  infinite  duration.  If  the 
duration  of  the  channel  realization  is  sufficiently  long,  one  can  think  to  coded-blocks  of  length  sufficient 
to  achieve  an  (almost)  error-free  transmission.  If  this  is  the  case,  it  is  known  that  the  probability  of  an 
erroneous  transmission,  for  a  given  transmission  rate  R,  is  dominated  by  the  outage  probability  on  that 
link,  i.e.  the  probability  that  Ch,  which  is  now  a  random  variable,  is  less  that  the  desired  rate. 

Pout^T  {Ch  <  R}  .  (71) 

In  this  section  we  derive  an  expression  of  the  outage  probability  of  the  system  described  in  section 
VI-.l  with  special  emphasis  on  the  geometry  of  the  system,  i.e.  on  the  relative  distances  between  the 
source,  the  relays,  and  the  destination,  and  on  the  power  balance  between  the  S2R  phase  and  the  R2D 
phase.  We  then  compare  the  obtained  result  with  the  outage  probability  of  a  system  which  employs  direct 
transmission  between  the  source  and  the  destination. 

In  order  to  obtain  the  desired  result,  we  will  first  calculate  the  outage  probability  in  the  S2R  phase, 
Pout,S2R,  defined  as  the  probability  that  at  least  on  of  the  ur  independent  links  between  the  source  and 
the  relays  is  in  outage  with  respect  to  R.  The  first  result  on  Pout,S2R  will  be  general,  in  the  sense  that 
it  will  be  a  function  of  the  hr  (possibly  different)  distances,  di,...,dk,  between  source  and  relays. 
Subsequently,  we  shall  assume  that  the  relays  are  located  at  the  same  distance,  dsR,  from  the  source,  and 
also  that  the  distance  between  the  relays  and  the  destination  is  (approximately)  the  same,  thus  referring 
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to  a  situation  as  the  one  depicted  in  Fig.  15.  This  assumption  will  simplify  our  successive  analysis  of  the 
outage  prohahility  in  the  R2D  phase,  Pout,R2D,  and  of  the  end-to-end  link,  Pout, tot,  that  we  carry  out 
in  paragraphs  (VI-.4.a)  and  (VI-.4.h).  Thanks  to  this  assumption  we  can  single  out  the  effect  of  both  the 
relative  distance  between  the  source,  the  relays  and  the  destination,  and  of  the  power  balance  between 
the  two  phases. 

4)  Source  to  relays  outage  probability:  Let  us  consider  the  link  between  the  source  and  one,  say  the 
k-th,  of  the  relays.  The  instantaneous  capacity  of  this  SI  SO  link,  measured  in  bits  per  channel  use 
(bpcu),  is  given  by  the  term  inside  the  expectation  operator  in  (61),  i.e. 

(72) 

=  log  >  /c  =  1, . . . ,  njj  (73) 

where  Vs  is  the  transmission  power.  For  this  link,  the  outage  event  for  a  given  “target  information  rate” 
of  R  bpcu,  is  defined  as 

Cs2R,k  <  R-  (74) 

It  is  easy  to  see  that  the  outage  events  (74)  can  be  expressed  as 

2 

4  <(2^-1)^,  k  =  l,...,nR.  (75) 

We  define  the  outage  probability  in  the  S2R  phase,  as  the  probability  of  the  union  of  the  outage  events 
on  the  single  links.  It  is  sufficient  that  one  of  the  links  is  in  outage  to  have  an  outage  in  this  phase. 
Equivalently,  an  outage  occurs  if  (at  least)  the  minimum  among  all  the  instantaneous  capacities  of  the 
S2R  links  is  less  than  the  required  rate  R,  i.e. 


Cs2R,k  —  log  1  + 


Vs  \hk\^ 

0-2  (l  -I-  d 


-Po«t,52itPl'{min(C'52it,l,  •  •  •  ,Cs2R,nii)  <  R}  ■  (76) 

Since  the  instantaneous  capacities  (72)  are  monotone  increasing  functions  of  the  respective  path  losses, 
we  can  equivalently  write 

Pout,S2R  =  Pr  I  min  (6, . . . ,  <  (2^  -  l)  ^  |  .  (77) 

Due  to  the  assumptions  on  the  (normalized)  fading  channel  coefficients  hk,  k  =  1, ,  ur,  we  have  that 
...  ,Ck  are  statistically  independent  exponential  random  variables  with  parameter 

Xk{l  +  dl)  ,  k  =  l,...,  UR.  (78) 

We  define  now 

^min(^i, . . .  ,^„^) .  (79) 
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It  is  known  that  the  minimum  in  a  set  of  independent  exponential  random  variables  is  itself  exponential, 
with  parameter  equal  to  the  sum  of  the  parameters.  Hence  ^  is  exponential  with  parameter 


Ur 


Ur 


A^Afc  -  ^  (1  +1 


(80) 


k=l 


k=l 


From  (77)  we  have 


Pout,S2R  =  Pr  <j  ^  <  (2^  -  l)  — 


=  1  -  e  (81) 

Let  us  now  assume  that  all  the  relays  are  located  at  the  same  distance  dsR  from  the  source,  i.e. 

dk  =  dsR,  fc  =  1, . . . ,  Ur.  (82) 


In  this  case  it  is  easy  to  show  that  the  outage  prohahility  (81)  becomes 

Pout,S2R  =  1  —  e 


^-nR{l+d^SR){2^-l)^ 


(83) 


In  the  sequel  we  will  use  this  expression  to  evaluate  the  end-to-end  outage  probability.  We  now  turn  our 
attention  to  the  R2D  phase. 

a)  Relays  to  destination  outage  probability:  If  the  data  sent  from  the  source  have  been  perfectly 
decoded  by  the  all  the  ur  relays,  we  end  up  with  a  set  of  terminals  that  have  the  same  data  and  may 
then  cooperate  to  transmit  them  to  the  destination.  In  practice,  the  set  of  relays  can  be  viewed  as  an 
antenna  array,  and  the  channel  between  this  array  and  the  destination,  as  a  MIMO  channel.  As  for  the 
SISO  channel,  we  consider  the  capacity  given  by  a  single  channel  realization: 

Cr2d  =  log  det  (l^  +  ^  w)  ,  (84) 

V  nRa^  J 

with  W  defined  in  (64). 

Under  the  simplifying  assumption  that  all  the  relays  are  at  the  same  distance  dRD  from  the  destination, 
the  diagonal  matrix  D  becomes  a  multiple  of  the  identity  matrix  of  size  hr,  i.e.  D  =  {l  +  '  Irr- 

In  this  case,  the  entries  of  the  matrix  HD  identically  distributed,  and  W  acquires  the  expression  in 
(67). 

It  is  hard  to  obtain  a  closed  form  expression  for  the  probability  density  function  of  the  instantaneous 
capacity.  However,  in  literature,  several  asymptotic  analyses  have  been  carried  out,  that  provide  a  useful 
mean  to  obtain  an  (albeit  approximate)  expression  of  the  outage  probability.  We  refer  to  the  asymptotic 
expression,  for  the  distribution  of  the  instantaneous  capacity,  derived  in  [59]  (see  also  the  reference 
therein). 
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More  specifically,  introducing  the  following  quantities 


L0\ 


la'^  (1  +  4d) 


Vr 


nR 


qo 


ro 


/3  -  1  -  ^2  +  ^ {P  -  1  - 


2uj 


1  —  /3  —  a;2  + 


2io 


(85) 

(86) 

(87) 

(88) 


it  can  be  shown  that  as  the  number  of  transmitting  and  receiving  antennas  goes  to  infinity,  and  the  ratio 
P  =  no/uR  tends  to  a  constant,  the  instantaneous  capacity  tends  to  be  distributed  as  a  Gaussian  random 
variable, 

Cr2d  ^  {i^c,o-q)  ,  (89) 


with  the  mean  and  variance  given  by 


fic  -  nR  {{1  +  P)  log  to  +  qoro  log  e 

+  log  ro  +  P log  (qo/P)}  (90) 

4-logelog(l-^).  (91) 

In  practice,  the  asymptotic  mean  and  variance  of  Cr2d  yield  a  close  approximation  to  the  statistics  of 
Cr2D  even  for  very  small  ur  and  riR. 

Therefore,  the  outage  probability  in  the  R2D  section  of  the  link,  with  respect  to  the  required  rate  R, 
can  be  approximated  with  the  Gaussian  Q-function,  as  follows: 

Pout,R2D  —  Q  ^  •  (^^1 

It  is  shown  in  [59]  that  this  expression  gives  a  close  approximation  to  the  actual  Pout,R2D  even  at  very 
small  values.  In  the  following  paragraph  we  exploit  the  results  on  the  two  phases  of  the  transmission  to 
obtain  an  expression  for  the  outage  probability  of  the  whole  system. 

b)  End  to  end  outage  probability:  The  end-to-end  communication  link  is  the  cascade  of  the  S2R 
and  the  R2D  links.  With  the  setup  of  section  VI-.l,  to  guarantee  that  the  information,  that  flows  through 
the  two  links,  can  be  perfectly  recovered  at  the  destination,  it  is  necessary  that  none  of  the  two  links  is 
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in  outage,  otherwise  the  information  is  destroyed  and  cannot  he  recovered.  Since,  for  a  given  rate,  the 
outage  events  in  the  two  links  are  independent,  it  is  easy  to  show  that  the  outage  prohahility  is  given  hy: 

Pout, TOT  —  Pout,S2R  P  Pout,R2D  Pout,S2RPout,R2D  •  (93) 

Assuming  that  the  relays  are  located  at  a  distance  dsR  from  the  source  and  on  the  line  connecting  the 
source  with  the  destination,  as  in  Fig.  15,  the  distance  between  the  relays  and  the  destination  is  given  hy 

dRD  =  dsD  -  dsR.  (94) 


Combining  (93),  (83),  and  (92),  we  obtain: 


Pout,TOT  =1  —  e 

+  Q 


nc  -  R 


C 


-  1  -e 


T  'Pat, 
“  'PS2R  cr'^ 


Q 


He  -  R 


(t: 


C 


(95) 

(96) 

(97) 


This  expression  depends  on  the  distances,  on  the  power  allocation,  represented  by  the  parameter  a,  and 
on  the  ratio  between  the  total  transmission  interval  and  the  sub-intervals  reserved  for  the  S2R  and  the 
R2D  phases.  Notice  that  now  /rc'  and  are  related  to  this  quantities  through  Vr  =  (1  —  a) 
and  dRD  =  dsn  —  dsR  ■ 


Fig.  17.  Outage  probability  of  the  cooperative  system  (colored  lines)  versus  the  outage  probability  of  a  direct  SIMO  link. 
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In  Fig.  17  we  have  plotted  the  outage  prohahility  obtained  with  the  relayed  transmission,  and  the  outage 
prohahility  of  the  direct  transmission.  It  was  assumed  that  the  destination  was  equipped  with  an  array  of 
n£)  =  3  antennas,  and  ur  =  6  relays  were  employed  for  cooperation.  The  distance  between  the  source 
and  the  destination  was  cIsd  =  100,  and  the  path  loss  exponent  t]  =  A.  The  ratio  between  the  average 
transmit  power  and  the  noise  variance  at  the  receivers,  'Pav/o'^,  was  100  dB,  and  the  time-sharing  factor 
was  Ts2r/T  =  0.9.  Different  curves  refer  to  different  values  of  the  power  allocation  parameter  a.  We 
can  see  that,  for  a  given  power  allocation  a,  there  is  an  optimal  location  of  the  relays,  whereas  for  a 
given  location  of  the  relays,  the  power  allocation  can  be  adjusted  to  obtain  the  best  gain  over  the  direct 
transmission  system. 

In  this  case  we  can  see  that,  for  a  relay  location  up  to  halfway  between  the  source  and  the  destination, 
there  exist  a  power  allocation  that  yields  a  significant  gain  over  the  direct  system. 

VII.  Conclusion 

In  conclusion,  distributed  space-time  coding  can  be  an  important  tool  to  reduce  the  overall  radiated 
power  in  wireless  networks.  We  have  considered  here  only  the  case  with  two  hops,  but  further  improve¬ 
ments  may  be  expected  in  the  multi-hop  case.  The  price  paid  for  these  advantages  is  the  additional 
signaling  required  to  coordinate  the  transmission  of  source  and  relay  nodes,  an  important  issue  which  is 
currently  under  investigation. 
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